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Introduced  as  a  corollary  of  the  strain  energy  density  theory  are  the  quanti¬ 
ties  dV/dA  and  dW/dV  which  represent,  respectively,  the  rate  of  change  of 
volume  with  surface  area  and  the  strain  energy  density  function.  ^Thoy 
getherj[  determine  the  energy  used  to  damage  a  differential  area-4A^in  the  pro¬ 
jectile  penetration  process.  ~Th^  orientations  of  the  damage  planes  form  the 
failure  path. 

A  numerical  procedure  is  developed  for  modeling  the  material  damage  pro¬ 
cess  during  projectile  penetration.  The  progressive  damage  pattern  for  each 
time  increment  is  exhibited  where  the  elements  fail  nonhomogeneously.  For 
blunt  projectiles  impacting  relatively  hard  targets,  the  ccyiditions  for  plug¬ 
ging  failure  are  met  soon  a f^er^  impact  with  ^v*ry>>  littreTtlow  of  material*  in---- ? 
the  radial  direction-  This*mode  w  failures  is  investigated  by  invoking  differ¬ 
ent  assumptions  in  the  state  of  the  failed  elements1: -ft ^he  present  model  can 
also  treat  the  phase  transformation  of  solid  where  shear  bands  are  formed  in 
regions  of  highly  localized  energy  states. 
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FOREWORD 


This  report  contains  the  complete  description  of  the  Axisymmetric  Dynamic 
Energy  Density  (ADED)  Code  which  was  developed  by  the  Institute  of  Fracture 
and  Solid  Mechanics  at  Lehigh  University  under  Contract  No.  DAAG46-83-K-0158 
supported  by  the  Army  Materials  and  Mechanics  Research  Center  (AMMRC).  Special 
acknowledgements  are  due  to  the  Project  Manager,  Mr.  J.  F.  Dignam  and  Technical 
Monitor,  Dr.  S.-C.  Chou  from  AW1RC.  Their  input  and  encouragement  have  been 
most  motivating  and  made  the  completion  of  this  work  possible.  The  authors 
also  wish  to  acknowledge  the  many  helpful  discussions  by  Dr.  J.  G.  Michopoulos 
of  the  Institute  of  Fracture  and  Solid  Mechanics  during  the  early  stage  of  this 
work. 

The  damage  concept  in  the  ADED  program  differs  fundamentally  from  all 
existing  codes  dealing  with  penetration  mechanics  in  that  the  theory  of 
plasticity  is  no  longer  employed  for  describing  the  permanent  deformation  of 
material.  Instead,  the  continuous  range  of  weakened  states  material  elements 
from  the  undamaged  to  the  completely  damaged  states  are  described  by  the  rates 
at  which  energy  is  dissipated  per  unit  area,  the  orientation  of  which  can  vary 
for  each  time  step  from  element  to  element.  The  full  range  of  uniaxial  strain 
rates  data  covering  the  impact  penetration  process  are  consistently  translated 
and  used  to  predict  multi -axial  stress  state  behavior  in  elements  that  had  at 
one  time  damaged  and/or  failed  when  subjected  to  compressive  stress  states  can 
also  be  distinguished. 

The  results  on  plug  formation  show  that  ADED  Code  is  fundamentally  sound 


and  can  predict  all  the  experimentally  observed  features  of  projectile 
penetration  damage  based  on  uniaxial  data  alone.  The  quantitative  accuracy 
depends  on  the  completeness  of  the  material  data  bank  and  the  finite  element 
mesh  size.  These  refinements  and  modifications  can  be  easily  incorporated 
into  the  ADED  Code.  There  is  no  doubt  that  the  program  can  be  further  stream¬ 
lined.  Nevertheless,  it  was  considered  essential  to  make  this  preliminary 
work  available  as  early  as  possible  to  those  who  have  had  a  long  standing 
interest  on  this  subject. 


DAMAGE  PREDICTION  OF  PROJECTILE  PENETRATION  PROCESS 
BASED  ON  ENERGY  DISSIPATION  RATE 


by 

G.  C.  Si h  and  D.  H.  Song 
Institute  of  Fracture  and  Solid  Mechanics 
Lehigh  University 
Bethlehem,  Pennsylvania  18015  USA 

ABSTRACT 

The  process  of  projectile  penetration  covers  a  wide  range  of  failure  modes 
depending  on  the  impact  velocity,  configuration  and  material  of  the  projectile 
and  target.  Such  behavior  has  not  been  adequately  described  by  the  classical 
continuum  mechanics  of  assuming  that  the  constitutive  relations  in  each  material 
element  are  known  as  a  priori.  This  difficulty  has  been  overcome  in  this  work 
by  application  of  a  new  concept  assuming  that  material  damage  occurs  nonhomo- 
geneously  throughout  the  target  and  can  be  uniquely  associated  with  the  rate  at 
which  energy  is  dissipated  in  a  unit  volume  of  material.  Introduced  as  a  corol¬ 
lary  of  the  strain  energy  density  theory  are  the  quantities  dV/dA  and  dW/dV  which 
represent,  respectively,  the  rate  of  change  of  volume  with  surface  area  and  the 
strain  energy  density  function.  They,  together,  determine  the  energy  used  to  dam¬ 
age  a  differential  area  dA  in  the  projectile  penetration  process.  The  orientations 
of  the  damage  planes  form  the  failure  path. 

A  numerical  procedure  is  developed  for  modeling  the  material  damage  process 
during  projectile  penetration.  The  progressive  damage  pattern  for  each  time  incre¬ 
ment  is  exhibited  where  the  elements  fail  nonhomogeneously.  For  blunt  projectiles 
impacting  relatively  hard  targets,  the  conditions  for  plugging  failure  are  met  soon 
after  impact  with  very  little  flow  of  material  in  the  radial  direction.  This  mode 


of  failure  is  investigated  by  invoking  different  assumptions  in  the  state  of  the 
failed  elements.  The  present  model  can  also  treat  the  phase  transformation  of 
solid  where  shear  bands  are  formed  in  regions  of  highly  localized  energy  states. 

INTRODUCTION 

Penetration  mechanics  has  been  a  subject  of  continuing  interest  because  of  its 
importance  in  military  application  for  developing  faster  projectiles  and  stronger 
armor.  The  early  works  in  this  area  are  mostly  empirical  that  involve  the  experi¬ 
mental  correlations  of  such  parameters  as  impact  velocity,  projectile  mass,  target 
thickness,  penetration  depth,  etc.  The  objective  has  bean  to  come  forth  with  some 
understanding  on  the  trade-off  between  the  pertinent  variables  that  govern  the 
scaling  of  models.  These  approaches,  however,  are  becoming  less  and  less  suitable 
when  applied  to  explain  modern-day  technology.  They  do  not  lei;d  themselves  to  any 
physical  insights  of  the  mechanics  of  penetration  and  can  be  costly  as  the  tests 
involve  too  many  variables.  The  advent  of  the  modern  computer  has  offered  many 
new  and  previously  untried  avenues  to  research  in  penetration  mechanics. 

During  the  past  two  decades,  the  U.S.  has  expended  considerable  efforts  toward 
the  formulation  of  sophisticated  computer  programs  [1-3]  to  explain  the  dynamic  re¬ 
sponse  of  materials  and  the  projectile-target  failure  phenomena.  There  is,  however, 
considerable  diversity  in  these  works,  particularly  in  the  application  of  failure 
criteria  and  constitutive  relations  that  involve  a  priori  assumption  on  material 
behavior  and/or  location  of  failure  path.  Lacking  in  particular  are 

(?)  a  unique  failure  criterion  that  can  consistently  explain  the  complete  ma¬ 
terial  damage  process  involving  the  transformation  of  solid  to  liquid  and/or  gas, 


(2)  thz  tKansZation  o 6  mzaAuAzd  matzfUaZ  pKopeAttzA  fa°m  4-un pie.  tzAtA  to 
muZtiaxAjdt  aV teA4  AtatzA  ok  nonhcmoge.ne.ouA  zneAgy  AtatzA  -in  thz  pznztKation  pkoczaa. 

The  majority  of  the  present-day  computer  codes  in  the  U.S.  have  failed  to  comply 
with  the  necessary  requirements  just  stated  and  hence  are  limited  in  their  predictive 
capability.  They  are  usually  developed  to  reproduce  the  experimentally  observed 
phenomenon  of  projectile  motion  and/or  target  failure  and  involve  many  empirical 
parameters  that  are  problem-specific.  This  is  indicative  of  the  fundamental  diffi¬ 
culties  associated  with  the  application  of  continuum  mechanics  for  explaining  fail¬ 
ure  that  involves  a  wide  range  of  energy  dissipation  rates  that  occur  in  the  pro¬ 
jectile  penetration  process.  Some  of  these  shortcomings  are  discussed  in  [4]. 

♦ 

Phenomenological  investigations  of  the  so-called  "adiabatic  shear  bands"  have 

been  observed  [5]  in  explosively  fragmented  shells,  impacted  plates,  projectiles, 

etc.,  and  studied  metal lographical ly  in  [6].  The  highly  localized  shear  strain 

5  7-1 

rate  was  estimated  to  be  of  the  order  of  10  to  10  sec  in  steel  such  that  phase 
transformation  of  the  metal  can  occur.  More  recently,  the  same  phenomenon  has  been 
observed  in  highly  but  slowly  compressed  metal  bar  specimens  [7-9]  where  ample  time 
was  available  for  heat  transfer  to  take  place.  White  shear  bands,  however,  still 
appeared  in  regions  of  localized  deformation.  It  is,  therefore,  important  to  dis¬ 
tinguish  the  transfer  of  heat  at  the  local  from  that  at  the  global  scale  level.  A 
quantitative  analysis  of  the  white  shear  band  problem  has  been  carried  out  in  [10] 
by  accounting  for  the  rate  at  which  energy  is  dissipated  to  damage  the  material 
during  loading.  It  was  shown  that  the  excess  energy  in  addition  to  yield  and  frac- 

The  term  adiabatic  was  invoked  to  describe  the  process  of  impact  occurring  so 
quickly  that  little  or  no  time  is  left  for  heat  transfer  to  take  place  between  the 
system  and  its  surrounding. 


ture  of  the  4340  steel  cylinder  in  compression  contributes  to  intense  local  heat¬ 
ing  that  can  lead  to  phase  transformation. 

The  main  objective  of  the  present  work  is  to  provide  a  methodology  that  can 
analyze  the  complete  projectile  impact  damage  process.  Incorporated  into  the  nu¬ 
merical  scheme  is  the  strain  energy  density  criterion  [4,11-13]  that  is  valid  to 
general  loading  conditions,  material  types  and  structure  configurations.  The  ini¬ 
tial  effort  is  to  develop  the  algorithm  for  carrying  an  analytical  description  of 
the  damage  process  leaving  out  the  details  of  shear  band  formation  and  material  in 
the  hydrodynamic  state.  To  reiterate,  the  basic  approach  is  able  to  describe  all 
failure  modes  depending  on  the  energy  dissipation  rates.  This  includes  the  trans¬ 
formation  of  solid  to  liquid  and/or  gas. 

COROLLARY  OF  STRAIN  ENERGY  DENSITY  THEORY 

The  selection  of  failure  criterion  for  describing  impact  damage  has  been  prob¬ 
lematic.  A  common  procedure  is  to  compare  experimental  data  with  theoretical  pre¬ 
dictions  based  on  different  assumptions.  This  approach,  however,  is  not  adequate 
because  the  differences  between  the  results  cannot  be  clearly  identified  with  phys¬ 
ics.  The  merits  of  any  failure  criterion  should  be  judged  by  its  versatility  and 

usefulness  in  explaining  a  wide  range  of  physical  phenomena  and  is  free  from  self- 
★ 

contradictions.  There  is  also  the  fundamental  problem  of  translating  the  nonlinear 
uniaxial  data  to  elements  in  a  structure  whose  stress  or  energy  states  vary  from 
one  location  to  another.  At  present,  the  only  widely  used  theory  for  describing 

★ 

The  maximum  normal  stress  criterion,  for  example,  contradicts  itself  when  applied 
to  the  running  crack  problem  where  the  maximum  stress  component  acts  parallel  to 
the  crack  plane  rather  than  normal  to  it  as  required  in  the  original  assumption. 


nonlinear  material  behavior  is  that  of  plasticity  where  yielding  is  considered 
to  be  the  mechanism  of  progressive  damage  before  fracture.  Difficulties  in  scaling 
the  grid  patterns  also  arise  when  applying  numerical  methods  such  as  finite  dif¬ 
ferences  and  finite  elements.  The  continuum  mechanics  theories  are  developed  by 
assuming  that  the  rate  of  change  of  volume  with  surface  area  for  each  element  tends 
to  zero,  i.e.,  dV/dA  0.  Such  a  condition,  of  course,  cannot  be  physically  real¬ 
ized  in  the  numerical  analysis  and  significant  errors  can  result  in  situations  where 
dV/dA  undergoes  large  gradients  due  to  inappropriate  selection  of  mesh  size  distri¬ 
bution.  Even  more  significant  is  the  role  that  dV/dA  plays  in  uniaxial  testing 
which  will  be  discussed  subsequently. 

Length  o£  Homogeneity.  In  the  mechanical  testing  of  uniaxial  specimens,  it  is 
necessary  to  define  the  resolution  of  measurements  in  terms  of  at  least  a  length 
parameter,  say  i  in  Figure  1(a),  that  describes  the  degree  of  uniformity  or  homo¬ 
geneity  of  the  stress  or  energy  state.  For  a  linear  response  of  the  uniaxial  stress 
and  strain,  the  parameter  dV/dA  remains  constant  while  dV/dA  changes  in  the  case  of 
nonlinear  behavior.  The  rate  of  change  of  dV/dA  with  time  is  indicative  of  pro¬ 
gressive  material  damage.  For  a  given  material  or  microstructure,  loading  rate 
and/or  specimen  size  may  be  altered  to  construct  a  data  bank  consisting  of  a  family 
of  curves.  Figure  1(b).  Each  point  on  these  curves  can  be  uniquely  identified  by 
defining  dV/dA  and  e  instead  of  o  and  e.  Since  dV/dA  is  a  geometrically  determina¬ 
ble  quantity  that  can  be  easily  computed  or  measured  for  an  uniaxial  specimen  and 
transferred  to  elements  in  a  complex  stress  state,  it  can  be  used  to  preserve  the 

it 

The  assumption  of  the  uniaxial  stress  and  strain  curve  to  coincide  with  the  ef¬ 
fective  stress  and  effective  strain  curve  is  inadequate  for  situations  where  dila¬ 
tation  also  contributes  to  failure  such  as  elements  near  the  crack  tip. 


plane  of  homogeneity  in  the  uniaxial  test  to  the  plane  on  which  damage  is  pre¬ 
dicted  to  take  place.  That  is  to  establish  a  relation  between  C(dV/dA)0,eQ]  and 
[(dV/dA) . ,e.J.  Let  a  be  the  angle  between  these  two  planes,  then 

COSa  (1) 

in  which  (dV/dA).  refers  to  homogeneity  associated  with  ith  plane  and  (dV/dA)Q  with 
a=0  coinciding  with  the  damage  plane  in  the  uniaxial  test.  There  remains  the  prob¬ 
lem  of  translating  the  uniaxial  strain  quantity  e  to  the  multiaxial  strain  state 
even  though  the  stress  quantity  a  is  no  longer  directly  involved  in  representing 
the  uniaxial  data.  This  involves  the  application  of  a  corollary  of  the  strain  en¬ 
ergy  density  theory  that  leads  directly  to  a  new  damage  theory  in  continuum  mechan¬ 
ics  that  includes  plasticity,  viscoplasticity,  etc.,  as  special  cases. 

A  Theoxy  06  Material  damage.  The  strain  energy  density  criterion  in  its  original 
form  [11-13]  assumes  that  the  strain  energy  density  function,  dW/dV,  varies  from 
one  location  to  another.  The  fluctuation  gives  rise  to  peaks  and  valleys  which  in 
mathematical  terms  are  the  stationary  values  of  dW/dV.  The  basic  postulate  is  that 
pKogxeAAive  matedal  damage  can  be  uniquely  aAi>o  dated  with  the  note  at  which  en¬ 
ergy  iA  dUAipated  in  a  unit  volume  ofi  mateKial.  In  other  words,  failure  modes  at 
different  scale  levels  are  assumed  to  be  uniquely  related  to  thresholds  of  dW/dV 
which  are  experimentally  measurable.  In  the  uniaxial  tests,  dW/dV  represents  the 
area  under  the  true  stress  and  true  strain  curve  depending  on  the  temperature  and 
moisture  level: 

if 

Homogeneity  can  change  from  element  to  element  in  a  solid  as  the  interaction  of 
material  with  loading  is  nonuniformly  distributed. 


37  =  /  aijdeij  +  f(AT.AC) 


with  AT  and  AC  being  changes  in  temperature  and  moisture  concentration  .  The 
stress  and  strain  components  are  denoted  by  o^.  and  e^.,  respectively. 

In  order  to  transfer  the  homogeneity  of  the  damage  plane  in  the  uniaxial  test 
to  the  prospective  plane  of  failure  of  an  element  in  multiaxial  stress  state,  it 
is  necessary  to  introduce  a  corollary  of  the  strain  energy  density  theory.  It  may 
be  stated  as  follows: 

The  hailwie  o £  an  element  it  at turned  to  occ.ua  on  the  ilk  plane  by  matching 
(dW/dA) .  mlh  meatuAable  uniaxial  data. 

For  an  element,  (dW/dA) .  can  be  written  as 


<$>,  ■  $>,<$> 


in  which  dW/dV  is  a  scalar  and  (dV/dA)^  is  proportional  to  the  slope  of  the  stress 
and  strain  curve.  Since  equation  (3)  must  be  unique  for  a  given  element,  the 
ith  direction  can  be  determined  from  the  condition 


v  **] 


m 


W  ■  const.,  i  =  ?,i 


In  two  dimensions,  equation  (4)  or 


Equation  (2)  shows  that  energy  can  be  stored  in  a  material  even  when  the  stresses 
are  zero.  Hence,  any  failure  criteria  based  on  stress  quantities  alone  are  neces¬ 
sarily  limited  in  application. 


r 1 

' J >  V 


v-\ 


/ dVx  _  ,dV> 

"  W 

5  n 


yields  the  angle  a  between  (s,n)  and  (x,y).  The  reference  damage  plane  is  as- 

* 

sumed  to  coincide  with  (x,y).  If  the  material  is  isotropic  and  homogeneous  , 
equation  (5)  reduces  to 


e_  =  e 
?  n 


These  are  not  the  principal  strains  but  those  acting  on  the  plane  with  the  same 
dV/dA  value  as  that  in  the  uniaxial  test.  In  view  of  equations  (5)  and  (6),  homo¬ 
geneity  of  the  multi  axial  stress  or  energy  state  can  be  described  by  [ (dV/dA) _  er] 

s »  s 

or  [(dV/dA)^  e^]  and  related  to  the  uniaxial  data  [(dV/dA)Q  eQ]. 

A  theory  of  material  damage  follows  immediately.  As  the  uniaxial  data  bank 
provides  known  values  of  (dV/dA).  and  dW/dV  for  each  energy  state,  the  actual  stress 
and  strain  path  of  each  element  can  be  derived  as  the  system  is  loaded  incrementally. 
The  procedure  for  constructing  the  stress  and  strain  history  of  a  typical  element 
is  illustrated  schematically  in  Figures  2  and  3. 

An  initial  stress  and  strain  of  a  given  element  is  assumed  for  the  first  load 

increment.  Equation  (6)  is  then  applied  to  find  a,,  (dV/dA),  and  e_  =  e  .  The 

I  ni 

predicted  damage  state  makes  an  angle  a-j  with  the  reference  state.  Figure  2(a). 

This  locates  the  point  "p,"  in  the  data  bank  with  coordinates  [(dV/dA)-., e  ]  as 

1  1  ?1 
shown  in  Figure  2(b).  The  next  increment  of  loading  gives  c^,  (dV/dA^  and 

=  e  ,  and  hence  the  point  "p9"  is  obtained.  Figures  3(a)  and  (b).  A  series  of 

*Expression  similar  to  equation  (6)  may  be  deduced  for  anisotropic  and  nonhomo- 
geneous  materials  in  which  case  the  relation  between  e  and  e  will  involve  mate¬ 
rial  constants. 


t  ■  ■  -- » '  A 
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(a)  Direction  of  homogeneity 
Figure  2.  Searching  of  point  p-j  in  data 


(b)  Location  of  point  p-j 
bank  for  the  first  load  increment. 


nl 


(b)  Location  of  point  p2 


Figure  3.  Searching  of  point  p2  in  data  bank  for  the  second  load  increment. 


points  Pj,  p£,  etc.,  are  found  until  the  complete  stress  and  strain  history  is 
derived.  Unloading  can  also  be  accomplished  incrementally  in  the  same  way.  Each 
element  will  follow  its  own  path  of  loading  or  unloading.  The  energy  dissipated 
in  damaging  each  element  can  be  assessed  with  accuracy. 

Preliminary  results  have  been  obtained  for  structural  members  undergoing 
static,  dynamic  and  fatigue  loadings.  They  indeed  exhibit  the  trend  anticipated 
on  physical  grounds.  In  the  case  of  a  slowly  growing  crack  in  a  stretched  plate 
[14],  the  equivalent  uniaxial  stress  and  strain  response  in  elements  near  the  crack 
tip  experiencing  more  dilatation  is  quite  different  from  those  away  from  the  crack 
where  the  plasticity  solution  gave  reasonable  results.  At  the  immediate  vicinity 
of  the  crack  tip,  large  deviations  are  expected  since  plasticity  accounts  only  for 
distortion.  Moreover,  the  damage  theory  predicts  a  much  sharper  rise  if  the  stress 
components  as  the  crack  tip  is  approached.  The  gradient  changes  for  each  increment 
of  crack  growth. 

FINITE  ELEMENT  FORMULATION 

The  finite  element  method  will  be  used  to  divide  the  projectile-target  system 
into  a  network  of  triangular  elements  such  that  each  element  possesses  three  nodes 
identified  by  a  convenient  coordinate  system.  In  the  case  of  axisymmetry  where  the 
projectile  has  a  circular  cross-section, the  (r,e,z)system  will  be  adopted.  Since  the 
result  does  not  vary  with  e,  it  suffices  to  consider  the  plane  r  and  z  as  shown  in 
Figure  4.  A  typical  triangular  element  with  nodes  i,  j  and  k  is  referred  to  in  a 
counterclockwise  direction.  A  composite  matrix  can  thus  be  formed  for  the  projec¬ 
tile-target  system  that  relates  the  displacements  of  the  nodal  points  of  each  ele¬ 
ment  to  the  external  forces  in  the  dynamic  structure.  Once  the  displacement  field 
is  known,  the  incremental  strain  can  then  be  evaluated  from  the  displacement  incre- 


Figure  4.  A  typical  triangular  element  in  rz-plane. 


ment.  At  this  point,  the  procedure  deviates  from  the  classical  approach  where  the 
stresses  are  obtained  from  the  strains  via  a  pre-assigned  constitutive  relation. 

The  newly  proposed  damage  model  utilizes  the  uniaxial  data  stored  in  the  bank  cover¬ 
ing  a  wide  range  of  strain  rates  that  may  differ  from  element  to  element  for  each 
increment  of  loading  or  projectile  advancement.  Based  on  the  damage  criterion  given 
by  equation  (3),  the  equivalent  uniaxial  stress  and  strain  history  for  each  element 
can  be  derived. 

ViAplcLceiMnti  and  S£ncU.nA.  For  problems  with  axial  symmetry,  the  displacement  vec¬ 
tor,  say  u,  depends  only  on  the  variables  r  and  z  as  shown  in  Figure  4.  Let  ur  and 
uz  be  the  displacement  components  in  the  r-  and  z-direction  of  u  associated  with  an 


element: 


The  displacement  vectors  of  the  nodal  points  i,  j,  k  denoted  by  u  ,  ir  and  u 
will  each  have  their  corresponding  components  in  the  r-  and  z-di recti on  as  fol- 


(u1}  >,  {uj> 


j  “H  k  “r 

J>  =<  >,  <uK>  =< 


For  the  triangular  element,  the  displacement  varies  linearly  with  r  and  z,  i.e.. 


{u(r,z)}  =  {c*i }  +  {cujlr  +  {a3>z 


in  which  {a^},  {agl  and  {a^}  are  column  vectors  that  can  be  expressed  in  terms  of 


{u1},  {uJ}  and  {u*},  i.e., 


(u(r,z) }  =  YU  J.  (ae+ber+cez){ue) 
e=i , j ,k 


such  that 


'i  ’  rjzk  -  Vj-  b(  *  zj  -  zk>  ci  -  rk  -  r) 


The  other  expressions  a-,  b.,  etc.,  can  be  obtained  by  cyclic  permutation  of  the 

J  J 

indices  i,  j  and  k.  The  quantity  D  stands  for  the  area  of  the  element  i,  j  and  k 
and  is  given  by 


L»J 


The  current  strain  e  at  time  t  will  be  calculated  from  the  incremental  strain 
Ae  at  At  as  follows: 

<;>t  ■  +  {a5}«  <13> 

where  is  the  strain  vector  at  time  t-At.  The  components  of  e  are 


(14) 


Starting  with  a  zero  state  of  initial  strain,  i.e.,  (e}Q  =  0,  the  strain  increment 
Ae  can  be  calculated  from  the  displacement  increment  Au  as  given  by 


» 

' a(aur) 

Aer 

3r 

Aez 

3(AUZ) 

3Z 

Ae0 

►  =  i 

AUr 

r 

AYr: 

r 

7 

3(AU  )  6(AU  ) 

+  4 

3Z  3r 

(15) 


Making  use  of  the  relations  in  equation  (10),  it  follows  that 


{Ae}  = 


I  be(AU  )( 

e«1,j,k  r 


_J.  ,  ce(Auz^ 

e-i > j  »k 


l  (^+  ba  +  ca  £)(auJ* 
e-itj.k  r  e  e  r;v  r' 


I  [c.(au  )e  +  b  (au  )e] 

e=1,j,k  e  r  e  2 


The  above  matrix  expression  may  be  written  as 


{Ae>  =  [Bi ,Bj,Bk] 


f  .  i 

{AU1} 


{AUJ’}1 


{AUk} 


provided  that  [Bg]  (e  =  i,j,k)  is  the  following  4x2  matrix: 


[Be]  = 


be  ,  0 

0  ,  ce 

ae  7 

—  +  b  +  c  —  ,  0 
r  e  e  r  * 

»  K 


Because  of  the  dependency  of  [Bg]  on  r  and  z,  the  strains  in  the 
longer  constant  as  in  the  problem  of  plane  extension.  Referring 
of  the  element  located  at 


i 

I 

k 

I 

b 

* 

i 

I 

(16) 


(17) 


(18) 


element  are  no 

IT  to  the  centroid 
e 
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equation  (17)  becomes 


|W}' 

(Ae)  =  [B.,  Bj,  Bk]|{AUj}  , 

{Auk} 

**  ~  . 

Equation  (20)  may  be  expressed  simply  as 
{Ae}  =  [B]{Aue> 

in  which 


(20) 


(21) 


[B]  *  [Bi,  Bj,  Bk]  (22) 

is  a  4><6  matrix.  With  reference  to  [1T3i  equation  (21)  yields  constant  strain  in 
the  element.  In  equation  (21),  ue  is  given  by 


S&vza&za  on  Damage  Plane.  The  nonhomogeneity  caused  by  the  variation  of  dV/dA 
throughout  the  projectile- target  system  can  strongly  affect  the  translation  of 
uniaxial  data  to  the  multi-axial  stress  state.  Equations  (5)  and  (6)  reveal  that 
the  normal  strain  components  control  this  effect.  Their  stress  and  strain  be¬ 
havior  must  be  derived  in  accordance  with  the  degree  of  nonlinearity  or  material 
damage  at  each  point  for  each  time  step.  More  specifically,  the  location  of  the 
damage  plane  determined  by  the  angle  a  in  the  rz-plane  must  be  found  from  equation 
(6)  for  a  material  that  is  initially  homogeneous  and  isotropic. 

The  stress  tensor  a  contains  only  four  components: 


{a}  =  l 


rz 


(24) 


They  correspond  to  the  strain  components  in  Figure  5(a).  Because  of  axial  sym¬ 
metry,  changes  in  the  9-direction  are  constrained  and  only  the  normal  strains  e 


and  ez  are  involved  in  the  adjustment  for  variation  in  dV/dA.  Referring  to  Figure 


5(b),  the  change  of  volume  with  respect  to  surface  area  in  the  £  and  n  direction 
* 


can  be  written  as 

(e„+£  )cOSa+£ 


(M)  - 
W  ■ 


9 


c,  ir 

_  (e  COSa+en )cOSa 
£  n  9 ' 


(25) 


and 


...  (e_+e  )c0Sa+eQ 

'(RT  „  "  ( e  .COSa+e. )COSa 

n  u  o 


(26) 


For  axisymmetry,  (dV/dA).  =  C(r,z)  and  eQ  in  equations  (25)  and  (26)  can  be  elimi 


nated  and  expressed  in  terms  of  C(r,z).  The  classical  plane  strain  condition  can 
be  realized  by  letting  Ol . 


VV,*y  r' 


Here,  (d V/dA)?  or  (dV/dA)^  is  related  to  (dV/dA)^  in  equation  (1)  through  cosa. 

As  a  result  of  the  corollary  of  the  strain  energy  density  theory,  equation  (5)  when 
applied  to  equations  (25)  and  (26)  leads  to  the  conclusion  given  in  equation  (6) 

or  the  condition  er  =  e  .  This  determines  the  orientation  of  the  damage  plane 

5  n 

through  a  such  that  the  appropriate  strain  rate  from  the  uniaxial  data  can  be  cor¬ 
rectly  transferred  to  each  element  in  the  projectile-target  system.  Note  that  the 
shear  strain  component  y^z  does  not  appear  in  equations  (25)  nor  (26).  It  is  there 
fore  related  to  xrz  in  the  usual  manner  through  the  shear  modulus  G. 

Once  [(dV/dA)^.,e^]  or  [(dV/dA)^,^]  are  determined  for  each  element  and  each 
time  step,  a and  can  be  found  from  the  material  data  bank  such  as  that  illus¬ 
trated  schematically  in  Figure  1(b).  Mathematically  speaking,  o or  can  be  ob¬ 
tained  from  the  relation 


a.  or  a  = 
5  n 


E(eS  or  en}’  eC  or  £n  -  ey 


e  or  e 

E(e?  or  en)  -  Bay[(  --)  -  1],  or  <  ey 


where  E,  ay  and  ey  are  respectively  the  Young's  modulus,  yield  strength  and  strain 
at  yield.  The  parameters  6  and  y  can  be  found  numerically  after  the  stress  and 
strain  for  a  given  element  has  been  traced  for  many  time  steps.  The  stress  compo¬ 
nent  a  is  related  to  eQ  by  an  expression  similar  to  that  given  by  equation  (27). 
The  stresses  in  the  triangular  elements  are  constant  as  a  consequence  of  constant 
strains  or  equation  (21). 


Nodal  FoticeA.  Once  the  stresses  within  each  element  are  found,  they  can  be  used  to 
determine  the  forces  acting  on  each  node  which  are  equilibrated  with  the  boundary 


tractions.  The  corresponding  nodal  displacements  can  then  be  obtained  by  lumping 
the  mass  of  the  element  at  the  nodes.  To  this  end,  the  principle  of  virtual  work 
will  be  employed.  Applied  will  be  the  condition  that  i&  the  nodes  are  -in  equi¬ 
librium,  the  total  virtual  work  done  by  the  forces  acting  through  any  arbitrary 
v-OUual  nodal  displacements  that  one  consistent  with  the  constraints  must  be  zero. 

Let  the  nodal  forces  be  denoted  in  matrix  form  as 


{ Fe }  =  < 


{F1} 

{FJ} 

{Fk} 


(28) 


For  an  arbitrary  nodal  virtual  displacement  Sue,  the  principle  of  virtual  work  may 
be  applied  to  yield 

T  -j-  j  t  n 

{6ue}  { Fe }  +  /  [{5e}  1  {cr}  -  {6u}'{f}]dV  -  /  {6u>'{T}dA  =  0  (29) 

V  -  '  ~  A  ~  ~ 

n 

in  which  f  and  T  are  respectively  the  body  force  and  traction  vector.  The  quantity 
<5e  is  the  virtual  strain  corresponding  to  the  virtual  displacement  <5u.  In  the  ab¬ 
sence  of  body  force  and  boundary  tractions,  equation  (29)  reduces  to 

{ 6Ue }T{ Fe }  =  -  /  {6e}T{(j}dV  (30) 

V  -  - 


With  the  aid  of  equation  (21),  it  follows  that 


{SU6}  {Fe}  =  -  /  {6ue}  [Bl'{0}dV 


(31) 


Since  { 6u  }  is  arbitrary  and  is  independent  of  the  volume  integral,  equation 
(30)  becomes 


{F6}  =  -  /  [I]T{a}dV 
v 


For  the  «.th  node,  the  force  vector  may  be  written  as 


(F*>  =  -  /  [B£]  {a} dV 
v 


in  which  [B is  given  by 


[b*]T  -  J* 


aQ+bQr+cQz 

be,  0  ,  * . f-*-.  ce 


0  ,  ce,  0  »  bft 


(32) 


(33) 


(34) 


The  quantities  D,  ag,  etc.,  are  the  same  as  those  defined  earlier.  The  volume  in¬ 
tegration  in  equation  (33)  may  be  carried  out  by  recalling  that  the  stresses  are 
constant  within  each  volume  element  and  applying  equation  (24)  the  nodal  forces 
components  F^,  F^,  etc.,  are  found.  With  2nrA  being  the  volume  of  an  element,  the 
expressions  for  F^,  and  F^  are 

Fr  -  -  »f[(Zj-zk)or  *  (rk-rj)trz]  -  jUo, 

(35) 

Fz  =  -  .f[(rk-rj)a2  *  (zrzk)V2] 

where  f2  =  0.  The  forces  on  the  other  nodes  such  as  F^,  F^,  etc.,  can  be  obtained 

9  r  z 

by  cyclic  permutation  of  the  indices. 


Equation  oi  Motion.  By  distributing  the  mass  of  the  element  evenly  at  the  three 
nodal  points  and  applying  the  nodal  forces  in  equation  (35),  the  acceleration  at 
the  nodes  may  be  obtained: 


^  Me 


(36) 


where  Me  stands  for  the  lumped  mass  at  the  eth  node.  The  velocity  vector  can  thus 
be  found  by  integrating  equation  (36).  Assuming  that  the  acceleration  is  constant 
for  a  small  time  increment  At,  then 


{Ue}  =  {0e>  +  {ue}.At 

'  t  '  t"  ~  Z 


(37) 


Referring  to  Figure  6,  (ue>  and  {ue }  are  the  velocities  just  after  and  before 

_1  t+  '  t‘ 

the  current  time  t  and  At  stands  for  the  average  time  increment  about  t,  i.e.. 


At 


(At)  +(At)  . 
t  t 
2 


(38) 


At  incipient  impact  t=0,  { ue }  are  the  initial  velocities  of  those  nodes  in 

'  t" 

the  projectile  which  come  in  contact  with  the  target  such  that  {ue>t  =  0  for 

(At)  .  Assuming  that  the  velocities  are  constant  over  the  time  increment  At, 
t" 

equation  (37)  can  be  integrated  to  yield  the  nodal  displacements 


{“ew  1  (“e>t +  <“e>t+<4t>t+ 


(39) 


where 


{ AU® }  =  {U®>  .(At)  , 

—  ~  X  '  X  * 


(40) 


is,  in  fact,  the  incremental  displacement.  The  nodal  coordinates  are  then  obtained 
from  the  displacements: 


(?  *t+At  =  *?*t  +  ^t+At 

To  summarize,  the  nodal  velocities  are  calculated  at  the  mid-point  of  the  time 
interval  which  represents  the  average  velocities  during  the  time  increment  At.  The 
nodal  displacements  and  accelerations  are  defined  at  the  beginning  or  the  end  of 
the  time  increment.  Refer  to  Figure  6  for  a  pictorial  representation  of  their  defi¬ 
nitions.  Use  is  made  of  the  "sliding-surface"  technique  developed  in  [15]  at  the 
contact  surface  of  projectile  and  target. 


(At) 


(At)  +- 
t 


time  t 


{“}t-At 


{U} 


{-}t+At 


- Displacement 


Figure  6.  Interpretation  of  displacement,  velocity  and  acceleration 
at  different  time. 
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Foc£ote  ConAide/uvtion.  By  application  of  the  conditions  given  in  equations  (5) 
and  (6),  the  orientation  of  the  damage  plane  defined  by  the  angle  a  in  Figure  5(b) 
is  first  obtained.  The  quantities  (dV/dA)^  or  (dV/dA)^  in  equations  (25)  or  (26) 
must  then  be  related  to  (dV/dA)Q  obtained  from  the  uniaxial  data,  say  in  Figures  1 
This  is  accomplished  by  application  of  equation  (1)  such  that  (dV/dA)..  in  equation 
(3)  stands  for  (dV/dA)^/cosa  or  (dV/dA)  ycosa.  At  a  particular  state,  say  p,  on 
the  true  stress  and  true  strain  curve  in  Figure  7,  dV/dA  being  proportional  to  the 
slope  can  be  expressed  by  a  length  parameter,  say  hp.  The  element  is  partially 
damaged  as  the  unloading  path  pq  will  not  coincide  with  the  loading  path  oy  which 
is  a  straight  line,  the  reversible  path.  The  shaded  area  oypq  is  the  energy  dissi 
pated  per  unit  volume  that  is  not  recoverable.  According  to  equation  (3),  the 
damage  is 


5  or  n 


.  u  /dWx 

■  hp  Wr 


such  that 


(42) 


fdWl 


<  M) 


<  h  rdwl 
_  -  hc  'dV' 

5  or  n  c 


(43) 


The  quantities  h  and  h  are  related  to 

y  c 

strain  curve  at  y  and  c  in  Figure  7  and 
areas  oyy1  and  oycc1,  respectively.  An 


the  slopes  of  the  true  stress  and  true 

(dW/dV )  and  (dW/dV)  correspond  to  the 

y  c 

element  is  assumed  to  have  failed  when 


fdwl 

W 


S  or  n 


h  A 
hc  V 


(44) 


This  condition  should  be  carefully  distinguished  from  that  in  equations  (42)  or 
(43).  In  the  present  analysis,  the  progressive  damage  of  all  elements  is  moni- 


True  Strain 

Figure  7.  True  stress  versus  true  strain. 

tored  at  all  time.  Those  reaching  the  condition  in  equation  (44)  are  considered 
to  be  failed  and  the  elements  will  be  removed  completely  from  the  finite  element 
grid  pattern  with  appropriate  adjustment  made  on  the  stress  state  due  to  changes 
in  geometric  configuration  and/or  the  initiation  of  cavities  or  cracks  in 
directions  dictated  by  the  strain  energy  density  theory. 

For  the  projectile-target  problem,  the  failed  elements  can  still  be  entrapped 

★ 

in  the  system  and  interact  with  the  others  as  the  damage  process  continues  .  In 
the  absence  of  strain  rate  data  on  the  projectile  and  target  material,  two  sets  of 
true  stress  and  true  strain  will  be  constructed  analytically  in  accordance  with 
the  generally  accepted  trade-off  relationship  between  yield  strength  and  fracture 
toughness.  These  two  cases  will  be  referred  to  as  Model  I  and  II. 

The  solid  may  transform  to  liquid  and/or  gas  at  high  impact  velocity  such  that  the 
rates  of  energy  dissipation  can  no  longer  be  adequately  described  by  the  area  under 
the  true  stress  and  true  strain  curve.  The  Hugonir;,  relation  expressing  the  hydro- 
dynamic  effects  of  material  may  come  into  play.  This  will  be  discussed  subsequently 


Model  I  -  The  projectile  and  target  material  are  assumed  to  possess  a  rela¬ 
tively  low  {ko-cAwkl  tougkneAA  for  the  range  of  strain  rates  considered. 

Model  II  -  The  projectile  and  target  material  are  assumed  to  have  a  rela¬ 
tively  high  ^fuictwiz  toughwAA  for  the  range  of  strain  rates  considered. 

The  corresponding  yield  stresses  or  stresses  at  which  permanent  deformation  oc¬ 
curs  are  kept  about  the  same  in  both  models.  The  time  leading  to  plugging  failure 
will  be  predicted  and  compared  so  that  the  sensitivity  of  experimental  strain 
rate  data  on  progressive  failure  can  be  demonstrated  and  better  understood. 

COMPUTER  ALGORITHM 

The  computer  algorithm  consists  of  three  main  portions.  They  are  referred  to 
as  INITI,  SEDOM  and  SPLOT.  Table  1  shows  the  flow  chart  for  INITI  that  deals  with 
the  initial  input  data  on  material  properties,  grid  generation  for  the  projectile- 
target  system  and  the  relevant  parameters  that  are  required  for  carrying  out  the 
failure  analysis.  The  chart  in  Table  2  gives  an  overall  view  of  SEDDM  that  incor¬ 
porates  the  strain  energy  density  theory  for  evaluating  the  damage  and  failure  of 
all  elements.  Sliding  nodes  are  introduced  at  the  interface  where  projectile  comes 
into  contact  with  the  target  in  order  to  ensure  displacement  compatibility  and  to 
avoid  overlapping  or  interpenetration  of  material  points.  The  damage  and  failure 
pattern  for  each  time  increment  can  be  exhibited  graphically  via  the  SPLOT  routine 
shown  in  Table  3.  Contours  of  constant  a  can  also  be  plotted  to  determine  the  path 
of  plugging  for  each  time  increment  At  which  is  chosen  according  to  the  interaction 
of  stress  waves  with  the  finite  element.  If  6  ,.n  denotes  the  minimum  dimension  of 
the  triangular  element,  then  At  should  be  smaller  than  the  time  required  for  the 

stress  wave  to  travel  across  6.  ,  i.e., 

min 


Table  2.  Flow  Chart  for  SEDDM:  Failure  Analysis 


The  accuracy  of  the  predictions  can  be  improved  by  making  the  element  size 
smaller  and  by  increasing  the  total  number  of  elements.  This  adds  complexity  to 
the  solution  as  the  computer  program  for  analyzing  the  progressive  damage  of  the 
projectile  and  target  already  involves  many  steps  where  the  stresses,  displacements, 
and  energy  densities  must  be  repeated  many  times  for  each  element  and  each  time 
step.  Refer  to  Tables  4,  5  and  6  for  a  more  detailed  account  of  the  computing  pro¬ 
cedures  involved  in  INITI,  SEDDM  and  SPLOT.  The  description  of  the  deck  cards  can 
be  found  in  [16]. 

IMPACT  OF  CYLINDRICAL  PROJECTILE 
ON  PLATE  TARGET:  MODEL  I  AND  II 

The  impact  of  a  cylindrically-shaped  projectile  on  a  finite  thickness  plate 
will  be  analyzed  by  application  of  the  axisymmetric  dynamic  code  that  incorporates 
the  strain  energy  density  material  damage  criterion  as  presented  earlier.  This 
program  shall  be  henceforth  referred  to  as  the  "Axisymmetric  Dynamic  Energy  Density 
(ADED)  Code".  The  impact  velocity  will  be  increased  incrementally  for  a  metal 
projectile  hitting  a  metal  target  until  perforation  occurs.  A  salient  feature  of 
ADED  is  that  material  properties  covering  strain  rates  from  10'^sec’1  to  lO^sec”1 
are  provided  in  the  computer  data  bank. 

This  covers  the  full  range  of  energy  dissipation  rates  that  produce  failure 
by  permanent  deformation,  spallation  and/or  fracture  in  the  form  of  plugging.  No 
a  priori  assumptions  are  made  on  the  mode  of  failure.  Unlike  all  the  other  codes* 

_ 

The  most  serious  limitations  in  all  these  codes  result  from  the  application  of 
the  von  Mises  yield  criterion  and  the  assumption  that  the  uniaxial  stress  and 
strain  data  coincide  with  the  effective  stress  and  effective  strain  used  in  the 
theory  of  plasticity.  Moreover,  they  cannot  consistently  account  for  the 
variations  in  strain  rate  effects  from  element  to  element. 


Table  4.  Block  Diagram  for  Computing  Procedures  of  INITI 


Table  5.  Block  Diagram  for  Computing  Procedure  of  SEDDM 
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Table  6.  Block  Diagram  for  Computing  Procedure  of  SPLOT 
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[1-3],  ADED  possesses  the  additional  ability  to  simulate  relative  damage  of  all 
elements.  This  allows  the  projectile  and  target  material  to  have  a  continuous 
range  of  weakened  states  from  the  undamaged  to  the  completely  failed  state.  Refer 
to  [16]  for  a  description  of  the  ADED  Code  where  the  card  deck  is  explained  in  de¬ 
tail  . 

PAjOjtatctz-TaAgzt  SyAtem.  Consider  the  projectile-target  system  illustrated  in 
Figure  8  where  both  the  projectile  and  target  are  made  of  4340  steel  with  a  Rockwell 
Hardness  number  of  52.  Referring  to  Figure  8,  the  geometry  is  such  that  Lp/Dp  * 

2.0,  ht  =  Dp  and  D^  =  10Dp.  These  proportions  coincide  with  the  example  problem 
in  [9],  The  finite  element  grid  pattern  is  shown  in  Figure  9  where  160  and  455 
elements  are  used  to  model  the  projectile  and  target,  respectively. 


Lp(44  units) 


D  (22  units) 

r 


Initial 

Impact 

Velocity  3  2,500  fps. 


Target  initially  at  rest 


D^. (110  units) 


Figure  8.  Schematic  of  projectile-target 

System:  Axially  symmetric  impact 


h t ( 22  units) 
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Figure  9.  Finite  element  grid  pattern: 

One-Half  Symmetry. 

MateAial  P/iope/ittea .  The  material  data  bank  in  the  ADED  Code  consists  of  eleven 
(11)  nonlinear  stress  and  strain  curves  modeled  by  the  relation 


a  *  Ee  -  80  [(f-)  -  1] 

y  y 


with  b  *  0.7879,  y  =  1.053  and  E  =  30  x  106psi.  The  yield  strength  ay  and 
corresponding  final  strain  are  given  in  Figure  10.  The  curves  labelled 
1,2,—  ,10  correspond,  respectively,  to  strain  rates  of  10  ,  10  ,---,10°sec"  . 

The  ADED  Code  traces  out  the  stress  and  strain  history  for  each  element  and  con¬ 
tinuous  damage  is  monitored  according  to  the  strain  energy  density  criterion, 
using  equations  (42)  and  (44).  Refer  to  Figure  11  for  the  relation  between  the 

yield  strength  a  and  critical  strain  energy  density  function  (dW/dV)  of  4340 
y  c 

steel . 


.•*V **.•’“  «*•«**  /•  .** 


“  *v*  “.**“> 


-35- 


Vp.  ^  ;v 


Model  I.  In  this  model,  the  failed  elements  are  considered  to  be  fractured  as 
they  reach  or  surpass  the  condition  in  equation  (44)  and  are  removed  completely 
from  the  analysis.  Figures  12  and  13  show  the  numbering  system  for  the  nodal 
points  and  elements.  The  results  will  be  described  in  a  series  of  time  steps 
labelled  1,2,  etc.  and  presented  only  for  those  time  steps  that  exhibit 
significant  damage  and/or  failure.  The  details  of  the  computer  program  are 
given  in  [16]  together  with  additional  numerical  results. 
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Figure  14  illustrates  the  projectile-target  damage  pattern  after  0.03053ysec 
of  impact.  The  corner  element  No.  160  on  the  projectile  has  already  failed  and 
been  removed  from  the  output.  Two  regions  are  identified.  The  dotted  elements 
nearest  to  the  contact  are  stressed  while  the  remaining  areas  are  unstressed  as 
the  waves  have  not  propagated  that  far.  At  t  =  0.04641ysec,  Figure  15,  the  two 
adjacent  rows  of  elements  with  numbers  151  ,  1 52 , —  ,159  and  201  ,  202,-209  that 
were  in  contact  have  failed  except  for  the  corner  elements  No.  210.  More  elements 
are  now  being  stressed  and  the  effect  of  wave  propagation  can  be  clearly  seen  in 
Figure  16  after  0.07715ysec.  The  elements  are  seen  to  fail  quickly  after  impact 
because  of  the  high  strain  rates  and  relatively  low  fracture  toughness  of  the 
target. 

The  values  of  a  locating  the  damage  plane  as  defined  in  Figure  5(b)  are  given 
in  Table  7  for  elements  near  the  contact  as  to  increases  from  0  to  0.4855ysec.  The 
direction  of  the  damage  plane  for  elements  No.  151  ,  152, —  ,160  at  initial 
contact  did  not  change  appreciably  up  to  failure.  For  elements  No.  210  and  211 
a  increased  slightly  i.e.,  the  axes  c  and  n  approach  towards  x  and  y  as  failure  is 


Time  =  0.04641ysec 


Table  7 


Angle  a  of  Damaged  Plane  for 
Element  Near  Contact  for  Model  I 


Element  No. 
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approached.  Oscillation  in  the  orientation  of  the  damage  plane  can  be  observed 
in  elements  No.  248  before  it  failed.  Table  8  gives  the  values  of  the  equiva¬ 
lent  uniaxial  stress  and  strain  in  the  e-  and  n-direction. 


Table  8.  Equivalent  Uniaxial  Stress  and  Strain 
in  Elements  Near  Contact  on  Damaged 
Plane  for  Model  I. 


Element  No. 

160 

151 

210 

Time 

(ysec) 

e=e  =e 
£  n 

(yin/in) 

a=a  =a 
£  n 
ksi 

e=e  =e 
£  n 
(yin/in) 

a=a=a 
£  n 
ksi 

£=e  =e 
£  n 

(yin/in) 

a-a  =a 
£  n 

ksi 

0.0050 

-2641 

-79.2 

-1710 

-57.3 

-784 

-23.5 

0.0105 

-5553 

-166.6 

-3596 

-107.9 

-1647 

-49.4 

0.0166 

-8756 

-262.7 

-5673 

-170.2 

-2596 

-77.9 

0.0232 

-12260 

-367.8 

-7955 

-238.7 

-3642 

-482.2 

-313.7 

-143.7 

0.0331 

(failed) 

-11360 

-340.8 

-5199 

-156.0 

0.0464 

(failed) 

-7145 

-214.3 

0.0611 

-271 .2 

(failed) 

The  corresponding  stress  and  strain  components  referred  to  the  r-  and  z-axis  can 
be  obtained  by  Mohr  circle  transformation.  Exhibited  graphically  in  Figure  17 
is  the  dynamic  stress  or  as  a  function  of  time.  The  normal  stresses 
in  element  No.  160  rise  more  sharply  with  time  and  it  fails  first.  Elements  No. 
151  and  152  in  the  projectile  and  elements  No.  210  and  211  in  the  target  failed 
subsequently.  The  corresponding  stresses  did  not  rise  as  sharply  with  time. 
Nonlinear  variation  is  not  observed  as  the  elements  failed  very  quickly 
on  account  of  low  fracture  toughness.  These  results  reflect  the  difference  in 
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Figure  17.  Equivalent  uniaxial  stress-time  curve  for 
elements  failed  at  contact  in  Model  I. 


Time  =  0.44352usec. 


□  Unstressed  Elements 
fl^  Stressed  Elements 


Figure  19.  Projectile-target  damage  pattern  after  0.44352ysec 
(Model  I). 
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3.742 
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2.563 
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3.168 
2.934 
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the  stress  and  strain  response  of  elements  owing  to  different  strain  rates.  This 
accounts  for  the  nonuniform  rate  of  energy  dissipation  from  element  to  element 
and  the  sequence  of  the  material  damage  and/or  failure  process. 

As  the  projectile  further  advances  up  to  t  =  0.09487ysec,  Figure  18  shows  that 
element  No.  211  being  under  element  No.  212  has  now  failed.  The  waves  have  also 
propagated  further  towards  the  back  side  of  the  target.  In  Figure  19  at 
t  =  0.44352ysec  the  interior  element  No.  286  is  seen  to  fail  as  a  result  of  dV/dA 
reaching  95.04  in  which  is  much  higher  than  the  three  neighboring  elements  with 
common  adjacent  sides.  They  are  elements  No.  249,  285  and  287  as  shown  in  Figure 
13  with  the  respective  dV/dA  values  of  2.384  in,  2.379  in  and  2.248  in.  Refer 
to  Table  9  for  some  typical  values  of  dV/dA  referring  to  elements  near  the 
contact  that  have  failed.  It  is  of  interest  to  note  that  the  value  of  dW/dV  = 
941.22  psi  in  the  unfailed  element  No.  249  is  higher  than  dW/dV  =  359.78  psi 
in  element  No.  286  which  has  failed.  This  is  due  to  the  high  strain  rate  effect 
as  dV/dA  is  proportional  to  the  slope  of  the  true  stress  and  true  strain  curve. 

The  dW/dA  criterion,  therefore,  includes  both  deformation  rate  through  dV/dA  and 
energy  dissipation  rate  via  dW/dV.  Hence,  failure  can  occur  at  locations  where 
dW/dV  may  be  low  but  dV/dA  can  be  high. 

Figure  20  plots  contours  of  the  damage  plane  orientation  with  constant  a. 

The  failed  element  No.  286  in  Figure  20  corresponds  to  contour  No.  5  with 

a  =  05.726°.  This  implies  that  the  5  and  n  axes  almost  coincide  with  the 

x  and  y  axes  along  which  uniaxial  data  are  taken.  It  is  essential  to 

recognize  that  the  direction  of  the  damage  plane  within  an  element 

★ 

need  not  coincide  exactly  with  the  path  of  fracture  or  plug  formation.  That 
★ 

It  is  common  experience  that  the  path  of  a  macrocrack  is  usually  assumed  to  follow 
the  general  formation  of  the  randomly  oriented  microcracks.  The  zig-zag  details 
are  microscopic  in  scale.  Their  influence  on  the  change  in  direction  of  the 
macrocrack  is  small  and  can  be  neglected  in  any  macroscopic  analysis. 


Time  =  0.48546ysec 


Figure  21.  Projectile-target  damage  pattern  after  0.48546ysec 
(Model  I). 


Orientation  of  Damage  Plane  at 
t  =  0.48546ysec. 


is  the  loca.i  iouiLuit  plant  and  global  ^nacXufit  path  do  not  necetiatuly 

coincide..  It  is  the  loci  of  the  failed  elements  that  form  the  path  of  fracture. 

Since  the  shape  of  the  triangular  element  is  chosen  arbitrarily,  it  cannot  be 

used  for  the  shape  of  the  cavity  or  crack.  The  condition  as  stated  in  equation 

★ 

(44)  corresponds  to  failure  say  permanent  deformation  prior  to  fracture.  A 
small  cavity  or  crack  is  thus  assumed  to  be  formed  at  element  No.  286.  As 
mentioned  earlier,  the  orientation  of  this  small  crack  does  not  have  to  coincide 
with  the  path  of  fracture  which  is  assumed  to  follow  the  contour  with  constant 
a  that  is  oriented  almost  vertically  through  elements  No.  285  and  322  and  side¬ 
ways  in  the  direction  of  elements  No.  248  and  249.  The  stress  state  in  these 
neighboring  elements  are  therefore  intensified  by  approximately  twenty  (20)  times 
in  accordance  with  the  well  known  factor  2a/b  where  a/b  is  the  geometric  aspect 
ratio  of  a  narrow  elliptical  cavity  or  crack.  Once  an  internal  cavity  or  crack 
is  formed,  fracture  continues  to  occur  quickly.  As  stress  waves  continue  to 
propagate  in  the  projectile  and  target.  Figure  21  shows  that  the  elements 
No.  248,  249  and  285  are  failed  at  t  =  0.48546ysec.  The  fracture  path  protrudes 
the  front  side  of  the  target  plate.  The  projectile  has  momentarily  lost  contact 
with  the  target.  Since  damage  and  failure  is  constantly  changing  process, 
the  orientation  of  the  damage  plane  for  each  element  also  changes  accordingly. 

★ 

The  common  notion  of  yielding  is  no  longer  applicable  in  the  newly  proposed 
damage  model  although  the  quantity  yield  strength  is  still  being  used.  It  is 
more  appropriate  to  refer  to  damage  by  permanent  deformation  as  the  concept  of 
yield  surface  is  not  needed. 

**  ... 

A  crack  is  defined  to  be  an  elongated  cavity  with  a  major  to  minor  axis  ratio 

of  ten,  say  2a/b  =  10.  The  longest  and  smallest  dimension  should  differ  at 

least  by  one  order  of  magnitude. 


Orientation  of  Damage  Plane  at 
t  =  0.56932usec. 
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Figure  26.  Constant  a  contours  after  0.56932usec 
(Model  I). 


The  contours  of  constant  a  in  Figure  22  are  seen  to  have  altered  as  compared 
with  those  shown  in  Figure  20  of  the  previous  time  step  with  t  =  0.44332psec. 

The  newly  failed  element  No.  285  is  now  aligned  along  contour  with  a  =  -5.66°. 
Without  further  stress  intensification,  element  No.  322  fails  on  the  next  time 
step  with  t  *  0.52734ysec,  Figure  23.  Fracture  is  seen  to  extend  directly  in 
the  negative  z-direction.  Shown  in  Figure  24  are  the  constant  a  contours  around 
the  fracture  path.  This  suggests  the  further  intensification  of  elements 
No.  321  and  358  extending  towards  the  back  side  of  the  target.  The  result  is 
shown  in  Figure  25  where  both  elements  No.  321  and  358  have  failed  together 
with  the  corner  element  No.  210.  The  corresponding  constant  a  contours  are 
shown  in  Figure  26.  Element  No.  357  failed  at  0.61136ysec  at  which  time  the  frac¬ 
ture  path  is  exactly  one-half  way  through  the  target  plate.  The  compressive 
wave  front  has  finally  reached  the  back  side  of  the  plate  and  reflection 
begins.  The  elements  with  damage  planes  of  a  *  -18.51°  is  almost  directly 
under  the  fracture  path,  Figure  28.  Element  394  remained  intensified  and 
failed  at  t  =  0.6114ysec.  As  the  target  continues  to  fracture,  the  stresses  in 
elements  No.  393  and  430  are  intensified  and  failed  when  t  =  0.69586psec.  This 
is  shown  in  Figure  29  together  with  the  constant  a  contours  in  Figure  30.  As 
the  plug  is  being  formed  elements  No.  429  and  466  are  further  stress  intensified. 
They  failed  at  t  =  0.73808usec  which  is  illustrated  in  Figure  31.  Figure  32 
gives  the  corresponding  contours  of  constant  a.  This  process  continues  by  fail¬ 
ing  the  next  two  elements  ,  No.  465  and  502  in  Figure  33.  The  constant  a 
contour  in  this  case  is  given  in  Figure  34.  The  failure  of  elements  No.  465 
and  502  occurred  at  t  =  0.94774ysec.  as  plugging  extends  further.  Figure  35. 

There  is  a  deviation  of  the  constant  a  contour  from  the  path  of  plugging.  This 
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Time  =  0.69586usec. 


”1  CZ1  Unstressed  Elements 
(y-jl  Stressed  Elements 


Figure  29.  Projectile-target  damage  pattern  after  0.69586ysec 
(Model  I). 
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Orientation  of  Damage  Plane  at 
t  =  Q.69585visec. 


(1)  38.80° 

(2)  29.43° 

(3)  17.58° 

(4)  5.72° 

(5)  -6.13° 

(6)  -17.99° 

(7)  -29.84° 

(8)  -39.21° 


indicates  the  tendency  of  spalling  as  the  waves  are  reflected  from  the  back 
side.  Refer  to  Figure  36.  The  finite  element  mesh  size,  however,  is  not 
sufficiently  refined  to  predict  the  details  of  the  local  failure  mode.  Complete 
perforation  is  completed  when  the  last  three  elemtns.  No.  501,  538  and  537 
failed  and  the  plug  is  dislocated  from  the  target  as  shown  in  Figure  37.  Orienta¬ 
tion  of  the  damage  planes  corresponding  to  complete  perforation  are  given  in 
Figure  38. 

In  this  example  of  a  blunt  projectile  impacting  a  relatively  hard  target 
made  of  4340  steel  with  a  Rockwell  hardness  number  of  52  and  low  fracture  tough¬ 
ness,  the  conditions  for  plugging  failure  are  met  very  soon  after  impact.  The 
flow  of  material  in  the  radial  direction  is  minimal  and  the  fracture  behavior  is 
very  brittle.  This  can  be  seen  from  the  sequence  of  diagrams  summarized  in 
Figures  39(a)  to  39(j)  inclusive  showing  the  initiation  and  completion  of  the 
plugging  process.  This  involved  nine  time  steps  from  t  =  0.4432ysec,  t  =  1.0723usec, 
that  cover  the  period  of  approximately  one-half  of  a  micro  second.  The  predicted 
events  involve  several  idealizations  that  can  be  easily  improved  in  future  cal¬ 
culations.  The  most  noticeable  inaccuracy  lies  in  the  size  of  the  triangular 
elements  that  resulted  in  unusually  wide  fracture  path  being  equal  to  one  tenth 
of  the  projectile  diameter.  This  also  led  to  the  overly  exaggerated  straight- 
line  shaped  plug.  A  refinement  of  the  element  meshes  along  the  prospective 
path  of  plugging  would  lead  to  a  more  realistic  prediction  of  plug  profile  and 
minimize  abrupt  change  in  fracture  path  as  indicated  by  the  missing  corner  in 
Figure  39(j).  The  relatively  low  fracture  toughness  of  the  material  data  in 
Model  I  also  contributes  to  the  brittle-like  failure  pattern.  These  short¬ 
comings,  however,  are  considered  to  be  of  minor  importance  in  contrast  to 
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Orientation  of  Damage  Plane  at 
t  =  1 .07226jjsec. 


(1)  40.62° 

(2)  30.86° 

(3)  18.51° 

(4)  6.16° 

(5)  -6.19° 

(6)  -18.54° 

(7)  -30.89° 

(8)  -40.65° 


Constant  a  contours  after  1.07226ysec 
(Model  I). 


(g)  t  =  0.7381vi$ec 


(h)  t  =  0.8642ysec 


the  major  improvements  achieved  in  the  ADED  Code  which  can  consistently  incorpor¬ 
ate  uniaxial  strain  rate  data  to  material  elements  in  the  projectile  and  target 
without  making  ad  hoc  assumptions.  Moreover,  many  of  the  drawbacks  in  the 
classical  theory  of  plasticity  that  are  inherently  embedded  in  all  previous 
computer  codes  are  overcome  in  the  present  damage  model  based  on  the  theory  of 
strain  energy  density. 

Further  insight  into  the  process  of  penetration  mechanics  can  be  gained  by 

studying  the  stress  and  strain  time  history.  Referring  to  those  elements  near 

the  region  where  contacts  were  made  at  initial  impact.  Tables  10  to  13  summarizes 

the  numerical  data  for  a.  a  ,  a.  and  x  in  a  series  of  time  steps  from 

z  r  u  rz 

O.OOSOysec  to  0.4855usec.  The  data  are  also  displayed  graphically  in  Figures  40 
to  43  inclusive.  It  is  immediately  noticeable  that  the  axial  stress  component 
for  those  elements  referred  to  in  Figure  40  is  compressive  up  to  failure.  The 
time  sequence  of  failure  for  elements  160,  151,  210,  211  and  248  is  clearly 
indicated.  The  corner  element  No.  160  on  the  projectile  failed  first  followed  by 
element  No.  151  at  the  center.  Failure  of  elements  No.  210  and  211  in  the  target 
occurred  next  with  element  No.  248  lacking  far  behind  as  it  is  further  away  from 
the  contact  surface.  The  amplitudes  of  the  radial  stress  component  ar  in 
Figure  41  are  much  lower  than  those  in  the  axial  direction.  Both  elements  No.  211 
and  248  in  the  projectile  experienced  tension  before  failure  as  the  nearby 
material  is  heavily  compressed  in  the  direction  of  the  projectile  motion.  The 
circumferential  expansion  of  the  projectile  and  target  is  evidenced  by  the 
variations  of  the  stress  component  a.  with  time  as  displayed  in  Figure  42. 

Shown  in  Figure  43  are  the  time  histories  of  the  shear  stress  component  xrz  in 
elements  near  the  contact  surface. 
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■140.2 

-90.8 

-90.8 

-41.6 

-41.6 

- 

•294.8 

-190.9 

-190.8  • 

.  -87.4 

-87.4 

-0.0004 

464.9 

-301.1 

-300.8 

-138.0 

-138.0 

-0.0021 

651.2 

-422.1 

-421.4 

-193.8 

-193.8 

-0.0075 

854.2 

-755.4 

-719.9 

-255.4 

-255.5 

-0.0224 

ail ed) 

-786.6 

-750.9 

-277.3 

-277.5 

-0.0248 

(failed) 

-926.7 

-381 .0 

-381.3 

-0.0605 

(failed 

-481.2 

-741 .1 

-0.213 

-563.4 

-844.6 

-1 .450 

(failed) 

-847.5 

-4.151 

(failed) 

-8.103 

-14.26 

-24.00 

-38.93 

-60.40 

-88.390 

-120.15 

-149.31 
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Table  11.  Radial  Stress  af  in  ksi  for  Elements 
Near  Contact  (Model  I) 


•.5<v 

•■V-V. 


Element  No. 


0.0050 

-18.3 

-11.8 

-11.8 

0.0105 

-38.4 

-24.9 

-24.9 

0.0166 

-60.4 

-39.3 

-39.3 

0.0232 

-84.4 

-55.2 

-55.2 

4 

-: 

4 

-1 

8 

-l! 

7 

-2 

1 

-3 

0.0949 

(failed)  129.9 

0.1144 

(faile 

0.1358 

0.1594 

0.1853 

0.0364 


0.0050 


0.0305 


Table  12.  Circumferential  Stress  a.  in  ksi  for 
Elements  Near  Contact  (Model  I) 


0  0 

014  -0 

067  -0 

193  -0 

-0. 

447  -0 

530  -0 

0.019 
0.088 
0.254 


106 

0 

i28 

0 

181 

0 

87 

0 

!22 

0 

47 

0 

0.0949 

(failed)  1.743  2.401 

0.1144 

(failed)  3.651 

0.1358 

5.228 

0.1594 

7.154 

0.1853 

9.399 

0.2138 

11.84 

0.2452 

14.24 

0.2798 

16.28 

0.3177 

17.60 

0.3595 

18.00 

0.4015 

17.61 

0.4435 

16.76 

0.4855 

310.00 

(failed 

Time 

(ysec) 


0.0050 

0.0105 

0.0166 

0.0232 

0.0305 

0.0331 

0.0464 

0.0611 

0.0772 

0.0949 

0.1144 

0.1358 

0.1594 

0.1853 

0.2138 

0.2452 

0.2798 

0.3177 

0.3595 

0.4015 

0.4435 

0.4855 


Table  13.  Shear  Stress  in  t  ksi  for 

rz 


Elements  Near  Contact  (Model  I) 


Element  No. 


21.4 

0.0 

0.0 

21.3 

0 

45.0 

0.021 

0.004-. 

44.7 

0.043 

0.0493 

70.8 

0.099 

0.016 

70.4 

0.199 

0.2285 

99.0 

0.287 

0.048 

98.3 

0.579 

0.6630 

129.5 

1.223 

0.188 

128.7 

1.345 

1.541 

(failed) 

1.352 

0.208 

139.5 

1.596 

1.823 

(failed) 

0.230 

196.0 

3.787 

5.255 

(failed) 

260.2 

16.580 

12.00 

333.8 

27.310 

23.81 

36.950 

39.85 

(failed) 

56.35 

71.92 

84.44 

91.09 

88.72 
74.64 
47.94 

10.72 
-31.64 
-69.06 
-98.59 
-23.86 
(failed) 


m 


'  vvv.v 


sec 
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Tables  14  to  17  Inclusive  outlines  the  strain  rates  in  the  axial,  radial 

and  circumferential  direction.  First  of  all,  strain  rates  in  the  axial  direction 

can  vary  by  several  orders  of  magnitude  from  10®  sec”1  in  element  No.  160  to 

10  sec-1  in  element  No.  248.  This  wide  range  of  strain  rates  was  covered  by  the 

data  bank  in  Figure  10.  As  it  is  to  be  expected,  the  strain  rates  did  not  vary 

as  much  in  the  radial  and  circumferential  directions.  They  varied  mostly  from 

10'1  sec”1  to  104  sec”1  except  for  the  shear  strain  component  in  element  No.  160 
5  -1 

that  reached  10  sec  .  In  retrospect,  it  is  essential  to  have  a  theory  such  as 
the  one  proposed  in  this  work  that  can  adjust  for  the  changes  in  strain  rates 
from  element  to  element. 

Model  II.  Consider  now  the  same  4340  steel  with  high  fracture  toughness  at  the 
same  rates  from  10”4  sec”1  to  10®  sec”1  as  shown  in  Figure  44.  The  areas  under 
the  true  stress  and  true  strain  curves  have  increased  as  compared  to  those  in  Fig¬ 
ure  10  for  Model  I.  The  relation  between  the  yield  strength  and  critical  strain 

energy  density  is  given  in  Figure  45  such  that  high  o  corresponds  to  low  (dW/dV) 

y  c 

and  vice  versa.  This  trade-off  property  is  typical  of  many  of  the  engineering  metal 
alloys.  The  ADED  Code  will  be  employed  to  analyze  the  damage  and  failure  pattern 
of  the  projectile  and  target  by  holding  all  other  variables  constant  other  than 
the  uniaxial  data  as  mentioned  earlier. 

With  the  tougher  material,  more  time  is  needed  to  damage  and/or  fail  the  target. 
Larger  time  steps  will  be  taken.  Moreover,  more  energy  will  be  dissipated  by  perma¬ 
nent  deformation.  Without  going  into  the  details  of  the  results  for  the  initial 
time  steps.  Figure  46  displays  the  damage  pattern  at  t  =  0.2287usec.  The  corner 
element  No.  160  is  broken  and  several  of  the  elements  in  contact  have  been  damaged 
and  they  are  shaded.  The  permanently  indented  surface  of  the  target  at  contact  be¬ 
comes  more  noticeable  in  Figures  47  and  48  corresponding  to  t  =  0.3643ysec  and 


0.0050 

-10.56 

-6.84 

-6.84 

-3.14 

-3.14 

- 

0.0105 

-10.60 

-6.86 

-6.85  . 

-3.14 

-3.14 

1.97x10 

0.0166 

-10.60 

-6.86 

-6.85 

-3.15 

-3.15 

-2.65x10 

0.0232 

-10.55 

-6.85 

-6.83 

-3.16 

-3.16 

-1 .33x10 

0.0305 

-10.45 

-6.82 

-6.78 

-3.18 

-3.18 

-4.45x10 

0.0331 

(failed) 

-6.84 

-6.79 

-3.18 

-3.18 

-3.47x10 

0.0464 

(failed) 

-6.58 

-2.94 

-2.94 

-1.84x10 

0.0611 

(failed) 

-2.58 

-2.58 

-6.11x10 

0.0772 

-1 .92 

-1.92 

-1.77x10 

0.0949 

(failed) 

-1.32 

-3.90x10 

0.1144 

(failed) 

-6.48x10 

0.1358 

-1 .04x10 

0.1594 

-1.59x10 

0.1853 

-2.27x10 

0.2138 

-2.99x10 

0.2452 

-3.54x10 

0.2798 

-3.67x10 

0.3177 

-3.09x10 

0.3595 

-1 .78x10 

0.4015 

-1 .98x10 

0.4435 

9.40x10 

0.4855 

1 .90x10 

Table  15.  Radial  Strain  Rate  er  (xlO^sec  ^ ) 
for  Elements  Near  Contact  (Model  I) 


Element  No. 


(nsec) 

160 

151 

152 

210 

211 

248 

0.0050 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0105 

2.84 

-0.87 

-1 .43 

2.82 

-0.18 

-0.18 

0.0166 

9.39 

-2.88 

-4.71 

9.30 

-0.61 

-0.61 

0.0232 

20.66 

-6.33 

-10.35 

20.37 

-1.34 

-1 .34 

0.0305 

37.76 

-11.59 

-10.95 

37.07 

-2.46 

-2.46 

0.0331 

(failed) 

-11.42 

-18.09 

36.75 

-2.46 

-2.46 

0.0464 

(failed) 

2.25 

16.98 

6.36 

6.33 

0.0611 

(failed) 

-14.63 

16.89 

16.74 

0.0772 

-63.43 

-85.42 

937.3 

0.0949 

(failed) 

-140.30 

-141  .5 

0.1144 

(failed) 

-889.3 

0.1358 

0.1594 

0.1853 

0.2138 

0.2452 

0.2798 

0.3177 

0.3595 

0.4015 

0.4435 

0.4855 


-316.9 

251.4 

75.81 

115.3 
141 .0 

153.4 
153.2 
133.8 

78.57 

-19.94 

-146.7 

(failed) 
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Table  16.  Circumferential  Strain  Rate  e.  (xlCrsec  ) 

u 

for  Elements  Near  Contact  (Model  I) 


WV 


Element  No. 

r.  •• 

r  V 

2  210  211 

248 

0.0050 

0.0 

0.0 

0.0 

0.0 

0.0 

- 

0.0105 

0.03 

-0.87 

-1.15 

0.03 

0.4 

0.20 

0.0166 

0.10 

-2.88 

-3.80 

0.90 

1.20 

0.66 

0.0232 

0.21 

-6.33 

-8.34 

1 .97 

2.64 

1 .46 

0.0305 

0.38 

-11.59 

-15.27 

3.61 

4.84 

2.69 

0.0331 

(failed) 

-11.42 

-15.05 

3.57 

4.78 

2.65 

0.0464 

(failed) 

-3.85 

2.29 

5.65 

5.10 

0.0611 

(failed) 

0.28 

6.79 

8.82 

0.0772 

1.49 

8.59 

14.77 

0.0949 

(failed) 

10.97 

18.34 

0.1144 

(failed) 

21.38 

0.1358 

24.51 

0.1594 

27.23 

0.1853 

28.85 

0.2138 

28.51 

0.2452 

25.51 

0.2798 

19.63 

0.3177 

11.60 

0.3595 

3.22 

0.4015 

-3.06 

Table  17.  Shear  Strain  Rate  erz  (xlO^sec  ^ ) 
for  Elements  Near  Contact  (Model  I) 


0.0050 

0.0105 

0.0166 

0.0232 

0.0305 

0.0331 

0.0464 

0.0611 

0.0772 

0.0949 

0.1144 

0.1358 

0.1594 

0.1853 

0.2138 

0.2452 

0.2798 

0.3177 

0.3595 

0.4015 

0.4435 

0.4855 


160 

151 

152 

210 

211 

248 

185.4 

0.0 

0.0 

184.7 

0.0 

- 

185.7 

1 .68 

0.03 

184.5 

0.34 

0.39 

185.2 

0.57 

0.09 

183.7 

1.12 

1.28 

183.5 

1 .23 

0.20 

182.1 

2.47 

2.83 

Vv\- 

180.3 

2.25 

0.38 

179.4 

4.54 

5.20 

(failed) 

2.23 

0.37 

179.6 

4.48 

5.12 

«*’  V 

(failed) 

0.16 

184.1 

7.13 

11.17 

(failed) 

190.0 

10.71 

19.97 

■‘vj 

198.1 

14.22 

23.81 

(failed) 

18.48 

39.23 

(failed) 

36.69 

31.47 

23.00 

11.11 

-3.60 

-19.44 

-33.52 

-42.47 

-43.94 

-38.60 

-30.46 

-21.38 

(failed) 
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Stress  (ksi) 


Time:  0.4075ysec. 


□  Unstressed  Elements 
(^1  Stressed  Elements 
gS  Damaged  Elements 


0.4075usec.  Elements  No.  159  and  158  of  the  projectile  are  fractured.  This  phe¬ 
nomenon  has  been  observed  experimentally  for  targets  made  of  hard  material.  The 
damage  zone  is  seen  to  increase  as  the  waves  propagate.  Figure  49  shows  that  the 

remaining  elements  in  the  first  layer  of  the  projectile  are  all  broken  except  for 

★ 

element  No.  154  acting  as  a  debris  or  fragment  .  The  debris  disappears  in  Figure 
50  at  t  =  0.4953usec  on  account  of  the  idealized  assumption  in  this  preliminary 
work  to  neglect  the  influence  of  the  failed  elements.  As  the  damage  zone  increases 
with  the  propagation  of  waves,  plug  initiation  is  seen  to  occur  at  element  No.  282 
in  Figure  51  when  t  =  0.5345usec.  Referring  to  Figure  13,  the  three  elements  that 
have  common  sides  to  element  No.  282  are  elements  No.  245,  281  and  283.  At  failure, 
the  rate  of  change  of  volume  with  surface  area,  i.e.,  dV/dA  in  element  No.  282  be¬ 
came  exceedingly  high  with  a  value  of  105.83  in  in  comparison  with  2.508  in,  2.566  in 
and  2.559  in  for  elements  No.  245,  281  and  283.  The  strain  energy  density  function 
dW/dV  in  elements  No.  245  and  281  are,  respectively,  20,977  psi  and  11,268  psi  both 
of  which  are  higher  than  the  value  of  dW/dV  *  8,826  psi  in  the  failed  element  No. 

282.  This  exhibits  the  nature  of  failure  at  high  strain  rates  where  the  area  under 
the  stress  and  strain  curve  alone  cannot  fully  describe  the  event  of  dynamic  fail- 

irk 

ure  .  The  constant  a  contours  at  plug  initiation  are  displayed  in  Figure  52.  The 
two  curves  with  a  =  -18.82°  and  -31.49°  are  directed  almost  vertically  from  element 
No.  282.  This  is  indicative  of  the  direction  of  plugging  that  is  headed  normal  to 


★ 

In  fact,  each  element  will  fail  at  a  different  time  and  can  be  analyzed  as  such 
if  refinements  in  time  steps  are  made.  The  influence  of  the  fragments  can  also  be 
treated  without  difficulty.  This  will  be  discussed  in  more  detail  subsequently. 

kk 

In  contrast,  dV/dA  does  not  change  appreciably  in  creep  loading  and  damage  can  be 
more  accurately  monitored  by  dissipation  due  to  dW/dV  alone.  It  becomes  obvious 
that  there  prevails  many  situations  where  both  dV/dA  and  dW/dV  can  play  equally  im¬ 
portant  roles  as  in  regions  near  the  tip  of  a  moving  crack  with  moderate  speed. 
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Orientations  of  Damage  Plane 
t  =  0.5345usec. 


(1)  41.90° 

(2)  31.88° 

(3)  19.21° 

(4)  6.53° 

(5)  -  6.14° 

(6)  -18.82° 

(7)  -31.49° 

(8)  -41.51° 


Figure  52.  Constant  a  contours  after  0.5345ysec  (Model  II) 


the  target  surface.  It  should  be  reiterated  that  the  global  fracture  path  in  the 
damage  theory  is  assumed  to  follow  the  direction  of  constant  a  line  and  not  a  it¬ 
self  which  has  been  emphasized  in  the  discussion  of  results  for  Model  I. 

Two  distinct  features  of  the  results  for  Model  II  are  observed.  First,  plug 
initiated  directly  under  the  corner  of  the  projectile  which  is  closer  to  experi¬ 
mental  observation  when  compared  with  Model  I  in  Figure  19  that  occurred  further 
away  in  the  radial  direction.  Materials  near  contact  are  severely  damaged  before 
failure.  This  did  not  occur  in  Model  I  as  the  fracture  toughness  of  the  material 
were  too  low  that  led  to  predominantly  brittle  fracture.  The  detail  features  of 
plugging  in  Model  II  can  still  be  improved  by  reducing  the  finite  element  mesh  size 
and  the  inclusion  of  failed  elements  in  the  calculation.  Refer  to  the  results  in 
Tables  18  and  19  for  some  typical  values  of  a  and  dV/dA  in  those  elements  near  the 
contact  surface  for  Model  II.  Table  20  summarizes  the  data  of  the  equivalent  uni¬ 
axial  stress  and  strain  curve  for  elements  No.  160,  151  and  210  on  the  damaged 
plane,  the  direction  of  which  is  given  in  Table  18.  Unlike  the  results  in  Figure 
17  for  Model  I,  the  curves  in  Figure  53  are  all  nonlinear  as  energies  are  dissipated 
in  addition  to  those  to  cause  fracture  or  failure.  Unloading  has  also  occurred  in 
element  No.  210.  This  is  shown  in  Figure  53. 

The  failure  of  the  internal  element  No.  282  creates  a  small  crack  and  intensi¬ 
fies  the  stress  state  in  its  neighboring  elements  No.  245,  246,  281  and  318.  A 

★ 

magnification  factor  of  20  is  employed  as  explained  in  Model  I.  This  leads  to 

* 

The  procedure  will  be  automated  in  the  future  ADED  Code  where  the  nodes  of  the 
triangular  elements  can  be  shifted  to  simulate  a  1/r  singular  field  for  the  strain 
energy  density  function  dW/dV  where  r  is  a  radial  distance.  This  character  being 
independent  of  the  constitutive  relations  of  the  material  is,  in  general,  valid 
for  simulating  the  effect  of  a  localized  defect  with  its  major  axis  differing  by 
an  order  of  magnitude  in  relation  to  its  minor  axis.  This,  however,  will  require 
the  use  of  isoparametric  elements  instead  of  the  constant  stress  elements. 
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Table  20.  Equivalent  Uniaxial  Stress  and  Strain  in  Elements 
Near  Contact  on  Damaged  Plane  for  Model  II 


Element  No. 


(ysec) 


e=e  =e 
£  n 

(yin/in) 

Hi 

-  10,730 

-322.0 

-  22,220 

-491.0 

-  47,370 

-502.2 

-  77,950 

-514.2 

-108,900 

-528.7 

t-e.-t 

n 

£ 

(yin/i 

-  6,913 

-  14,460 

-  30,930 

-  48,840 

-  68,070 

-  78,490 

-  89,630 
-100,100 
-110,000 
-118,700 

(failed) 


-207.4 

-433.9 

-471.0 

-478.7 

-486.3 

-490.1 

-494.2 

-522.2 

-525.6 

-528.6 


e=e^=e 

(yin/i 


-  3,149 

-  6,611 
-13,320 
-15,920 
-13,940 
-10,260 

-  6,146 

-  2,930 

-  1,089 

-  1,614 


-  94.5 
-198.3 
-399.7 
-439.9 
-380.4 
-270.2 
-146.7 

-  50.2 

5.0 

-  10.7 


Time  (ysec) 


Figure  53.  Equivalent  uniaxial  stress  and  strain  for  elements 
near  contact  (Model  II). 


the  failure  pattern  in  Figure  54  at  t  =  0.5977ysec  where  fracture  is  extending  in 
both  the  upward  and  downward  direction.  Elements  No.  245,  246,  281  and  318  are 
now  all  failed.  The  waves  have  propagated  to  the  back  side  of  the  target  plate 
and  spread  radially.  Following  the  constant  a  contour  in  Figure  55,  elements  No. 
209,  210,  317  and  354  are  stress  intensified  while  plugging  continues.  Figure 
56  shows  the  resulting  fracture  path  where  plugging  has  extended  through  the 
front  surface  where  elements  No.  209  and  210  failed  at  t  =  0.6500ysec.  At  this 
time,  element  No.  208  is  still  intact  but  is  on  the  verge  of  being  failed.  The 
damage  zone  has  spread  to  the  mid-thickness  except  for  element  No.  347  as  indi¬ 
cated  in  Figure  56.  The  corresponding  constant  a  contours  are  given  in  Figure  57 
from  which  elements  No.  353  and  390  are  seen  to  be  situated  along  the  path  of 
prospective  plugging.  At  t  =  0.6900ysec,  plugging  has  penetrated  more  than  one 
half  of  the  target  thickness  as  elements  No.  317  and  354  fail.  A  cyl indrically- 
shaped  plug  is  formed  as  shown  in  Figure  58  where  elements  No.  208  and  381  are 
about  to  fail.  The  fracture  path  is  directed  along  the  constant  a  contour  No.  (7) 
in  Figure  59.  Elements  No.  353  and  390  are  thus  intensified.  As  the  fracture 
path  continues  to  extend  by  breaking  elements  No.  389  and  426  as  shown  in  Figure 
60  at  t  =  0.7340usec,  more  elements  are  damaged.  The  reorientation  of  the  damage 
plane  is  displayed  in  Figure  61.  Again,  contour  No.  (7)  with  a  =  -31.16°  is  di¬ 
rected  almost  normal  to  the  target  surface.  This  angle  did  not  differ  appreciably 
from  that  of  contour  No.  (7)  in  Figure  59.  Elements  No.  425  and  462  are  next  in¬ 
tensified  and  fractured  as  time  increases  to  0.8942ysec,  Figure  62.  Failure  also 
occurred  in  element  No.  242  at  the  top  of  the  plug.  The  constant  a  contours  ahead 
of  the  fracture  path  in  Figure  63  are  now  slightly  slanted  in  line  with  elements 
No.  461  and  498.  It  is  interesting  to  note  that  as  compressive  waves  are  reflected 
from  the  backside  elements  No.  453,  456,  457  and  460  leading  ahead  are  first  dam¬ 
aged  while  elements  No.  418,  419,  422  and  423  further  back  are  still  undamaged. 


Time:  0.5977ysec. 

Unstressed  Elements 
Stressed  Elements 
Damaged  Elements 


Figure  54.  Projectile-target  damage  pattern  after  0.5977ysec  (Model  II). 


-105- 


WWW 


Orientations  of  Damage  Plane 
t  *  0.6500ysec. 


Time:  0.6900usec. 


Figure  58.  Projectile-target  damage  pattern  after  0.6900ysec  (Model  II) 


Figure  59.  Constant  a  contour  after  0.6900usec  (Model  II) 
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Time:  0.7340visec. 


60.  Projectile-target  damage  pattern  after  0.7340usec  (Model  II) 
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Orientationsof  Damage  Plane 
t  =  0.8942ysec. 


(1)  39.43° 

(2)  29.70° 
17.39° 

5.08° 
■  7.23° 

(6)  -19.54° 

(7)  -31.85° 

(8)  -41.58° 


(3) 

(4) 

(5) 


\  \ 
'(7)  (6) 


(3) 

(4) 

(5) 


Figure  63.  Constant  a  contours  after  0.8942usec  (Model  II) 


Figure  67.  Constant  a  contours  after  1.1004ysec  (Model  II) 


Figure  64.  This  exhibits  the  tendency  of  spallation.  At  this  time,  the  majority 
of  the  constant  a  contours  in  Figure  65  are  aligned  almost  normal  to  the  plate 
surface.  The  final  stage  of  plugging  is  completed  when  elements  No.  497,  534  and 
533  fail  at  t  *  1.1004ysec.  Figure  66  shows  that  element  No.  463  next  to  the 
fracture  path  is  also  damaged.  The  plug  is  dislocated  from  the  target  and  is 
damaged  both  externally  and  internally.  The  failure  of  elements  No.  419  and  422 
inside  the  plug  suggests  fragmentation  as  the  plug  leaves  the  target  with  the 
exit  velocity.  Figure  67  summarizes  the  orientations  of  the  residual  damage  planes 
in  both  the  projectile,  target  and  plug. 

The  sequence  of  events  starting  from  plug  initiation  to  dislocation  is  shown 
in  eight  (8)  time  steps  as  outlined  in  Figures  68(a)  to  68(h).  For  the  tougher 
material  used  in  Model  II,  plug  initiation  started  at  a  later  time  with 
t  =  0.5345nsec  as  compared  with  t  =  0.4432ysec  for  Model  I  in  Figure  39(a).  More¬ 
over,  the  location  was  moved  to  an  element  directly  under  the  corner  of  the  pro¬ 
jectile  rather  than  at  an  outward  radial  distance.  There  were  also  significant 
differences  in  the  damage  and  failure  patterns  of  the  projectile- target  system 
prior  to  plug  initiation.  A  glance  at  Figures  68  reveals  that  the  plug  experienced 
more  damage  in  Model  II.  As  a  consequence,  more  time  was  also  required  for  full 
perforation.  Since  both  Model  I  and  II  contain  several  idealized  assumptions,  the 

results  must  be  interpreted  accordingly.  Aside  from  the  coarseness  of  the  finite 

★ 

element  mesh  size,  damping  force  were  not  included  in  the  present  analysis.  This 
done  purposely  in  order  not  to  cloud  the  findings  with  any  arbitrariness  and/or 
artificial  effects.  It  can  be  accounted  for  in  future  analyses  where  damping  de- 


•k 

This  is  usually  added  into  the  computer  scheme  by  introducing  the  so-called  "arti¬ 
ficial  viscosity",  a  quantity  that  is  not  known  and  assumed  arbitrarily. 
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pends  on  load  history  and  varies  nonhomogeneously  in  the  system.  The  projectile 
velocity  obviously  becomes  unrealistically  high  when  this  factor  is  ignored. 
Nevertheless,  what  has  been  achieved  for  the  first  time  is  an  analytical  demonstra¬ 
tion  of  how  uniaxial  data  affect  the  failure  modes  arising  from  projectile  penetra¬ 
tions. 

Summarized  in  Tables  21  to  24  are  the  time  dependent  stress  components  az,  ar, 
oa  and  t  _  for  elements  No.  151,  156,  160,  210,  211  and  248.  These  locations  are 
chosen  to  illustrate  the  different  degree  of  stress  levels  prior  to  failure.  The 
corner  element  No.  160  of  the  projectile  is  most  severely  loaded  at  initial  impact. 
This  is  seen  in  Figure  69  where  the  curve  corresponds  to  the  axial  stress  az  in 
element  No.  160  which  dropped  sharply  as  it  is  intensified  in  compression.  It  then 
levels  off  before  failure.  Elements  No.  151  and  156  that  also  come  into  contact 
with  the  target  behaved  in  a  similar  fashion  but  failed  at  a  later  time.  For  those 
elements  No.  210,  211  and  248  in  the  target,  the  axial  stress  az  attains  an  oscil¬ 
latory  character  as  they  are  damaged.  The  component  or  in  these  same  elements 
also  varied  appreciably  with  time  and  is  shown  in  Figure  70.  The  trend,  however. 
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Table  23.  Circumferential  Stress  aQ  in  ksi  for 
Elements  Near  Contact  (Model  II) 


Element  No. 


Time 

(usee) 


0.0200 
0.0420 
0.0928 
0.1543 
7 
6 


Time  (ysec) 


No.  248  \  No.  211 


210,  211  and  248  for  Model  II. 


Element  No.  210 


is  similar  co  az  except  for  the  curves  corresponding  to  elements  No.  151  and  156. 
They  oscillated  slightly  with  time.  As  a  result  of  axisymmetry,  the  change  in  the 
circumferential  direction  is  expected  to  be  relatively  small.  Figure  71  shows 
that  o  varied  only  slightly  with  time.  Similar  results  for  the  shear  stress  com- 
ponent  trz  is  given  in  Figure  72. 

The  nonuniform  strain  rates  at  the  different  locations  can  be  evidenced  from 

the  data  in  Tables  25  to  28  for  the  ez,  er,  e0  and  erz>  They  vary  from  10®  sec~^ 

2  -1 

to  10  sec  and  change  from  one  location  to  another.  It  is  clearly  seen  that  the 
strain  rate  variations  in  element  No.  160  are  the  largest  as  it  fails  first.  The 
changes  in  element  No.  248  that  are  further  away  from  the  initial  impact  are  not  as 
appreciable. 

N0NH0M0GENE0US  ENERGY  DENSITY  DISSIPATION  RATE 

It  cannot  be  overemphasized  that  the  failure  behavior  of  all  materials  are 
intimately  associated  with  the  rates  at  which  energy  is  dissipated  per  unit  volume 
or  unit  area.  For  problems  that  involve  a  wide  range  of  time  scale,  both  quanti¬ 
ties  dV/dA  and  dW/dV  can  be  equally  important.  Plug  initiation  was  found  to  occur 
at  locations  of  very  high  dV/dA  whereas  dW/dV  was  relatively  low  comparing  with 
those  in  the  neighboring  elements.  This  is  indicative  of  the  rate  effect  being 
proportional  to  the  slope  of  the  uniaxial  stress  and  strain  curve  that  is  related 
to  dV/dA.  When  the  loading  rates  are  very  low  such  as  those  experienced  in  creep 
deformation,  changes  in  dV/dA  are  small  and  the  rates  of  energy  dissipation  are 
governed  mostly  by  dW/dV.  In  general,  the  corollary  of  the  strain  energy  density 
theory  [4]  as  stated  by  equation  (3)  provides  a  more  consistent  treatment  since 
no  a  priori  assumption  is  made  on  the  form  or  nature  of  the  constitutive  relations 
This  is  why  conceptual  difficulties  often  arise  when  applying  the  dW/dV  criterion 


Table  26.  Radial  Strain  Rate  er  (x  10^  sec"^) 
for  Elements  Near  Contact  (Model  II) 


•v./v* 


mm 


Element  No. 


0.0200 
0 
8 
3 
7 
6 


160 

151 

210 

2 

0 

0 

0 

4.68 

-  1.40 

4.56 

-  0 

12.83 

-10.35 

12.75 

-  8 

32.14 

-17.06 

32.77 

-  4 

29.22 

-17.04 

37.40 

-27 

(failed) 

-14.98 

-  0.15 

-51 

-15.40 

-  49.01 

-71 

‘Mv.-: 
• , » » » 
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Table  28.  Shear  Strain  Rate  erz  (x  10"*  sec"^) 
for  Elements  Near  Contact  (Model  II) 


Element  No. 


Time 


isec) 


160 


151  210  248 


0200 

188.4 

0 

185.6 

- 

0420 

180.2 

2.75 

178.6 

6.2 

.0928 

176.9 

13.34 

166.2 

42.0 

1543 

214.4 

19.02 

181.9 

88.4 

2287 

280.8 

15.22 

179.5 

78.2 

2716 

(failed) 

7.78 

188.6 

44.6 

3188 

1.25 

186.4 

4.3 

3643 

-  3.38 

192.4 

-14.7 

in  conjunction  with  the  theory  of  plasticity.  The  yield  condition  of  von  Mises 
considers  only  the  distortional  component  of  dW/dV  while  the  relative  maximum  of 
dW/dV  that  is  assumed  to  coincide  with  locations  of  excessive  distortion  also  in¬ 
cludes  the  dilatational  component  of  dW/dV.  Such  a  conflict  will  always  prevail 

★ 

if  the  failure  or  damage  criterion  is  not  inherently  embedded  in  the  stress  or 
strain  analysis. 

To  illustrate  the  nonhomogeneous  character  of  dW/dV  and  dW/dA,  elements  No. 

15,  152,  160,  210,  211  and  248  are  selected  with  their  variations  as  a  function 
of  time  given  in  Tables  29  and  30  for  Model  I.  Because  of  the  brittle-like  fail¬ 
ure  behavior,  most  of  the  energies  are  dissipated  in  terms  of  failure  by  fracture. 
Displayed  in  Figures  73  and  74  are  respectively  the  time  dependent  character  of 
dW/dV  and  dW/dA.  The  curves  for  elements  No.  151,  156,  160,  210  and  211  all  rise 
very  sharply  up  to  the  point  of  failure  except  for  element  No.  248  that  follows  a 
much  more  gradual  change.  Since  both  the  variations  of  dW/dV  and  dW/dA  with  time 
are  very  similar,  the  influence  of  dV/dA  is  seen  to  be  small.  Significant  change 
in  dV/dA  occurred  only  at  plug  initiation. 

Tables  31  to  33  outline  the  results  on  energy  density  as  a  function  of  time  for 
Model  II.  Again,  the  difference  between  dW/dV  and  dW/dA  shown  in  Figures  75  and 
76  is  small  as  dV/dA  in  general  remained  fairly  constant  except  at  the  time  of  plug 
initiation.  The  curves  for  elements  No.  151,  156  and  160  in  Figures  75  and  76  in¬ 
creased  monotonical ly  with  time  while  those  for  elements  No.  211  and  248  achieved 
oscillation.  This  is  because  both  dW/dV  and  dW/dA  contain  both  elastic  and  dissi¬ 
pated  energy.  The  elastic  portion  can  increase  or  decrease  with  time  depending 

on  wave  propagation.  The  dissipated  portion  of  dW/dV  or  (dW/dV)p,  however,  must  al- 
★ 

Stress  analysis  is  a  misnomer  in  the  strain  energy  density  theory  because  the 
method  determines  the  stress  distribution  via  the  quantity  dV/dA  and  uniaxial  data. 
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Table  29.  Strain  Energy  Density  Function  dW/dV 
in  psi  for  Elements  Near  Contact 
(Model  I) 


Table  30. 


Strain  Energy  Per  Unit  Area  dW/dA  in 
lb/ in  for  Elements  Near  Contact 
(Model  I) 


Element  No 


152  to  156 


41 .710 


22.760 


(failed)  25,280 


0.1144 

(failed)  351.8 

0.1358 

580.5 

0.1594 

821.4 

0.1853 

999.4 

0.2138 

1,051.7 

0.2452 

1,010.6 

0.2798 

1  ,049.1 

0.3177 

1,460.6 

0.3595 

1,854.5 

0.4015 

2,224.9 

0.4435 

2.514.7 

(failed) 
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Table  31 


Strain  Energy  Density  Function  dW/dV 
in  psi  for  Elements  Near  Contact 
(Model  II) 


v:v.v?d 


S'V  *?  *» 


Element  No. 


160  151  156  210  211  248 


0 

0 

0 

0 

0 

- 

11.5 

0.4 

0.8 

0 

0 

55.7 

24.9 

23.6 

4.20 

2.108 

0 

111.9 

53.2 

48.7 

13.31 

3.195 

190.3 

84.0 

78.0 

20.62 

4.811 

0 

(failed) 

100.9 

94.1 

25.23 

4.811 

o.r 

for  elements 


ways  rise  monotonically.  This  is  shown  in  Figure  77  where  the  slope  of  the  curves 
gives  an  indication  of  the  severity  of  damage.  Element  No.  160  has  the  largest 
slope  in  the  (dW/dV)p  versus  time  plot  and  fails  first.  These  curves  yield  infor¬ 
mation  on  the  rate  at  which  energy  is  dissipated  as  a  function  of  time. 


REMARKS  ON  REFINEMENT  OF  ADED  CODE 


As  mentioned  throughout  this  report,  this  work  on  the  development  of  ADED 
is  considered  to  be  preliminary.  The  emphases  have  been  placed  on  applying  the 
basic  concept  of  the  strain  energy  density  theory  [4]  for  predicting  projectile 
penetration  damage  from  uniaxial  data  alone.  This  required  a  change  in  one  of 
the  fundamental  approaches  of  continuum  mechanics;  that  is,  the  constitutive  re¬ 
lations  for  the  material  elements  will  no  longer  be  assumed  as  known  but  will  be 
derived  for  each  element  and  time  step.  Such  a  feature  is  necessary  for  describing 
the  nonhomogeneous  rates  of  deformation  in  penetration  mechanics.  Although  there 
remains  a  number  of  refinements  that  should  be  incorporated  into  ADED,  they  are, 
however  manageable  and  not  regarded  as  overwhelming.  A  description  of  these  re¬ 
finements  will  be  given. 

The  most  serious  drawbacks  in  many  of  the  presently  developed  computer  codes 
on  penetration  mechanics  are  inherent  in  the  use  of  the  classical  theory  of  plas¬ 
ticity  together  with  arbitrarily  assumed  failure  criteria  that  can  lead  to  incon¬ 
sistencies  and  contradictions.  Compatibility  between  failure  and/or  damage  cri¬ 
terion  and  stress  analysis  is  necessary  for  making  reliable  predictions.  This 
has  been  accomplished  through  equations  (3),  (42)  and  (44)  by  assuming  that  the 
energy  dissipated  per  unit  area  in  a  given  time  interval  can  be  uniquely  identi¬ 
fied  with  the  state  of  material  damage.  The  predicted  results  on  plugging  for 
Model  I  and  II  are  obviously  idealized  because  of  the  simplifying  assumptions 
which  can  be  overcome  without  difficulties. 

Energy  Vluipcution  R ateA.  The  ADED  Code  can  easily  accommodate  higher  rates  of 
energy  dissipation  that  can  no  longer  be  adequately  described  by  the  area  under 


the  uniaxial  stress  and  strain  curves  via  dW/dV.  Information  on  pressure  P  ver¬ 
sus  volume  V  can  be  used  instead.  A  more  suitable  quantity  is  the  internal  en¬ 

ergy  density,  say  U.  The  Mie-Griineisen  equation  of  state  is  a  case  in  point: 

r 

P  2  (a-|U  +  &2u2  +  a3^3)0  -  2^0  +  rU(l+y)  (47) 

in  which  p  =  VQ/V  -  1,  a^  (i  =  1,2,3)  are  material -dependent  coefficients  and  r 

is  the  Grtineisen  coefficient. 

Predictions  on  the  appearance  of  shear  band  [10]  during  plugging  can  also  be 
made  by  generalizing  equation  (3)  into  the  form 

<w>T  1  <w>T  <w>T  '48> 

at  a  given  temperature.  A  much  more  refined  finite  element  mesh  is  needed  for 
this  purpose.  Equation  (48)  can  be  uniquely  identified  with  the  energy  dissipa¬ 
tion  rates  at  which  phase  transformation  of  the  metal  takes  place. 

Finite.  Element  Size.  The  size  of  the  finite  elements  in  the  prospective  region 
of  failure  and/or  damage  can  be  refined  so  that  more  realistic  shape  of  the  plug 
can  be  predicted.  This  can  be  easily  done  in  future  calculations. 

FnagmentcutionA .  Those  elements  that  have  failed  as  a  result  of  impact  between 
the  projectile  and  target  can  remain  as  individual  particles  in  the  subsequent 
calculations.  The  sliding  surface  technique  can  be  applied  such  that  each  of  the 
fragments  can  transmit  momentum  and  energy  while  displacement  and  stress  continuity 
will  no  longer  be  required.  In  this  way,  they  continue  to  contribute  in  the  dam¬ 


age  process. 


FKa.ctuA.z  TrUticLtLon.  As  internal  elements  are  fractured,  voids  or  small  cracks 

are  created  causing  local  stress  or  energy  intensification.  This  can  be  handled 

automatically  by  using  the  isoparametric  nodes.  The  mid-nodes  in  a  triangular 

element  may  be  shifted  to  one-quarter  of  the  distance  from  the  point  where  the 

1/r  character  is  desired  in  the  strain  energy  density  field  where  r  is  a  radial 

distance.  This  simulates  the  character  of  an  internal  flaw  regardless  of  the  ma- 
★ 

terial  behavior  .  This  invariant  character  of  the  strain  energy  density  function 
makes  the  failure  analysis  consistent  through  the  loading  history. 

Nonh.omogzne.ou6  Vamping.  The  predicted  projectile  velocities  are  unrealistically 
high  for  both  Model  I  and  II.  This  is  because  no  artificial  viscosity  has  been 
introduced  that  prevails  in  all  other  computer  codes.  It  is  felt  that  such  a 
factor  will  not  be  constant  and  should  not  be  incorporated  arbitrarily.  It  can 
vary  from  location  to  location  for  each  time  step  depending  on  the  rate  of  energy 
dissipation  which  is  highly  nonhomogenzous.  The  quantity  (dW/dV)p  in  equation 
(42)  can  be  calculated  to  define  an  effective  damping  force  Fp  that  depends  on 
the  space  variables  and  time.  This  effect  when  included  into  ADED  will  lower  the 
projectile  velocity  in  a  cumulative  fashion. 

Time  I ncn.eme.nt  ok  Step.  The  selection  of  time  increment  or  step  in  the  numerical 
computation  is  essential  and  is  intimately  associated  with  the  element  size  and 
location.  These  variables  must  be  compatible  with  the  rates  at  which  energy  is 
being  nonuniformly  dissipated  while  the  material  is  damaged  and/or  failed.  The 
fluctuation  of  the  strain  energy  density  function  yields  information  on  how  the 
time  and  space  variables  should  be  scaled. 

* 

The  singular  character  of  the  stress  depends  on  the  constitutive  relations  and 
would  not  be  a  suitable  choice  in  this  approach  where  the  stress  behavior  changes 
from  element  to  element  for  every  time  step. 


ThAe.t-VmuM>ioYit> .  It  is  obvious  that  the  (dW/dA)^  concept  expressed  by  equation 
(3)  can  be  easily  extended  to  problems  of  oblique  impact  where  no  symmetry  pre¬ 
vails.  The  determination  of  the  damage  plane  for  an  homogeneous  and  isotropic  ma¬ 
terial  would  then  involve  equating  (dV/dA)?,  (dV/dA)^  and  (dV/dA)^  in  the  three 
orthogonal  directions  £,  n  and  ?.  The  finite  element  grid  pattern  would  be  con¬ 
siderably  more  complicated  but  the  basic  procedure  remains  unchanged. 

Many  of  the  aforementioned  refinements  are  already  being  incorporated  into 
ADED  [16].  In  principle,  the  method  applies  equally  well  for  describing  the  ir¬ 
reversible  behavior  of  material  in  the  solid,  liquid  or  gaseous  state.  It  is  the 
versatility  and  generality  of  the  strain  energy  density  concept  that  makes  the 
theory  attractable. 
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